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ABSTRACT 


This  research  presents  a  systematic  approach, to  the  optimal 
design  of  spatial  structures  for  minimum  veight  subject  to  con¬ 
straints  on  stress  and  geometry.  The  optimization  procedures  dis¬ 
cussed  are  general  end  may  be  applied  to  structures  which  can  be 
analyzed  by  matrix  displacement  or  finite  element  methods. 

Two  methods  of  mathematical  programing  are  applied  to  obtain  a 
minimum  weight  design.  The  first  is  the  sequential  unconstrained 
minimization  technique  (SUMT),  and  the  second  is  the  method  of  con¬ 
strained  steepest  descent  with  state  equations  (CSDS),  Both  of  these 
techniques  require  derivatives  of  the  objective  and  constraint  func¬ 
tions  to  improve  estimates  of  the  optimum  design.  In  many  structural 
problems,  it  is  very  difficult  or  impossible  to  compute  these  deriv¬ 
atives  exactly;  existing  structural  analysis  algorithms  are  generally 
not  equipped  to  compute  these  derivatives.  In  order  to  take  full 
advantage  of  existing  analysis  capability,  the  programming  techniques 
in  this  research  have  been  developed  assuming  that  such  derivatives 
are  not  available. 

Optimal  structural  design  problems  are  characterized  by  an  objec¬ 
tive  function  (the  veight),  state  variables  (the  stresses  and  deflec¬ 
tions),  design  variables,  state  equations  (the  structural  analysis), 
and  constraints  which  may  be  fuactions  of  the  design  and  state  vari¬ 
ables,  When  the  state  equations  are  used  to  write  sill  of  the 
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constraints  as  functions  of  the  design  variables,  a  nonlinear  pro¬ 
gramming  problem  results.  The  sequential  unconstrained  minimization 
technique  reduces  the  constrained  nonlinear  programming  problem  t-o  a 
sequence  of  unconstrained  problems  which  can  be  solved  using  existing 
unconstrained  minimization  techniques.  A  SUMT  program  was  written 
for  this  research  using  Powell's  method  of  unconstrained  minimization 
without  derivatives.  The  required  minimization  of  a  function  along  a 
line  uses  a  combination  of  a  Fibonacci  search  (to  bracket  the  minimum) 
and  a  quadratic  approximation  of  the  minimum. 

The  method  of  constrained  steepest  descent  differs  from  the  usual 
nonlinear  programing  problem  in  that  the  state  equations  and  the 
state  variable  constraints  appear  explicitly  in  the  formulation.  This 
provides  a  natural  matching  of  the  essential  features  of  the  design 
problem  and  the  method  used  to  obtain  its  solution.  The  design  prob¬ 
lem  is  linearized  about  a  candidate  design  and  the  desired  improvement 

in  the  design  variables,  6x,  is  required  to  be  small  by  demanding  that 
T  -1  ■> 

6x  w  6x  =  £*,  where  £  is  a  small  number  and  w  is  a  positive  definite 
weighting  matrix.  The  Kuhn-Tucker  necessary  conditions  are  then 
applied  to  the  resulting  nonlinear  problem.  As  a  direct  consequence, 

6x  is  specified  in  terms  of  two  components;  5x^  which  reduces  the 
objective  function  consistent  with  the  constraints,  and  6Xg  which 
directs  the  search  for  a  minimum  back  to  the  feasible  region  if  con¬ 
straints  have  been  violated.  The  method  was  applied  using  both  exact 
and  approximate  derivatives,  so  that  its  effectiveness  when  derivatives 


1 

; 

i 

? 

i 


are  not  available  could  be  assessed 


A  sp&tidl  stimcture  which  occurs  frequently  in  practice  is  the 
plane  frame  with  out-of -plane  loads.  Although  such  structures  are 
generally  made  up  of  relatively  few  members,  they  may  have  many 
design  variables  since  several  design  parameters  must  be  specified 
for  each  member.  The  programing  methods  were  applied  to  a  number 
of  two  at>d  three  member  frames  of  this  type.  Fro®  the  results,  it 
appears  that  CSDS  has  significant  advantages  over  SUMT  both  in  terms 
of  computational  time  and  the  number  of  times  that  candidate  designs 
must  be  analyzed.  The  results  also  show  that  CSDS  performs  as  well 
when  derivatives  are  approximated  S3  it  does  when  they  can  be  computed 
exactly.  The  effectiveness  of  SUMT  is  reduced  significantly  if  the 
derivatives  are  unavailable. 

Abstract  approved;  .  t  dissertation  supervisor 

*  fyab  ■  title  department 
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Key  Words  —  Optimum  design,  minimum  weight,  structures,  spatial 
structures,  mathematical  programming,  spatial  frames. 
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CHAPTER  I 


INTRODUCTION 


1.1  The  Class  of  Problems  to  be  Solved 

Hie  subject  of  optimum  structural  design  has  received  much  atten¬ 
tion  in  recent  years.  The  optimal  design  of  simple  structural  elements 
and  structures  consisting  of  these  elements  has  been  studied  exten¬ 
sively.  A  number  of  investigators  have  considered  plane  frames  and 
trusses,  but  spatial  frames  have  received  little  attention.  It  is  to 
this  topic  that  this  research  is  directed.  The  problem  to  be  solved 
is  the  minimum  veight  design  of  a  spatial  frame  subject  to  constraints 
on  stress  and  geometry. 

One  type  of  spatial  structure  that  occurs  frequently  in  practice 
is  the  plane  frame  with  out-of-plane  loads.  Similar  frames  are  often 
found  in  automotive,  construction  vehicle,  and  agricultural  equipment 
applications.  Real  design  problems  of  this  type  initiated  this 
research  and  are  used  as  example  problems  for  the  solution  techniques 
investigated.  Although  such  structures  are  generally  made  up  of  rela¬ 
tively  few  members,  they  may  have  many  design  variables  since  several 
design  parameters  must  be  specified  for  each  member.  In  addition  these 
frames  are  often  required  to  support  or  transmit  loads  at  many  points 
in  the  structure.  When  such  frames  are  mass  produced  the  design  which 
requires  the  least  material  has  a  significant  economic  advantage . 
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The  design  must  also  satisfy  constraints  on  stresses  and  geometry.  In 

moat  cases  the  structure  will  fail  by  fatigue;  therefore,  the  maximum 

stresses  in  the  frame  must  be  restricted  to  values  well  below  the 

elastic  limit  of  the  material.  Limits  on  the  design  variables  of  the 

frame  may  be  dictated  by  current  manufacturing  capabilities  or  other 

i 

related  limitations.  It  is  apparent  from  this  discussion  that  struc¬ 
tural  design  problems,  like  almost  all  design  problems,  have  a  merit 
oi  objective  function,  a  set  of  design  variables,  a  set  of  state  vari¬ 
ables,  such  as  stresses  or  deflections  which  describe  the  behavior  of 
the  frame,  a  set  of  state  equations  which  determine  the  state  at  a 
given  design  point,  constraints  on  the  design  variables,  and  con¬ 
straints  on  the  state  variables. 

Two  methods  of  mathematical  programming  are  applied  to  obtain  a 
minimum  weight  design.  The  first  is  the  sequential  unconstrained 
minimisation  technique,  and  the  second  is  a  constrained  steepest 
descent  method  which  uses  the  state  equations  directly  in  the  optimi¬ 
sation  process.  Both  of  these  techniques  require  derivatives  of 
objective  and  constraint  functions  to  predict  better  approximations  to 
the  optimum  design.  These  derivatives  may  be  ciaiberscoe  or  Impossible 
to  compute  exactly.  In  addition  it  is  desirable  in  structural  optimi¬ 
sation  to  take  advantage  of  existing  analysis  algorithms  which  in  gen¬ 
eral  yield  only  function  values,  not  derivatives.  In  order  to  take 
full  advantage  of  existing  analysis  capability,  the  programming  tech¬ 
niques  should  be  effective  when  only  function  values  are  available. 

The  optimisation  procedures  discussed  in  this  research  are  general 
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and  nay  be  applied  to  spatial  fr ernes  analyzed  by  matrix  displacement 
or  finite  element  methods. 

In  Chapter  II  the  method  of  structural  analysis  is  discussed  and 
the  formulation  of  the  optimization  problem  as  a  mathematical  pro¬ 
gramming  problem  is  presented.  The  sequential  unconstrained  minimi¬ 
zation  technique  is  outlined  in  Chapter  III  and  the  constrained  steep¬ 
est  descent  method  with  state  equations  is  developed  in  Chapter  IV. 

In  Chapter  V  the  application  of  the  programming  methods  to  seme  example 
problems  is  discussed  and  the  results  for  these  examples  ivre  presented. 
Conclusions  concerning  the  results  and  the  relative  merit  of  the  design 
methods  are  discussed  in  Chapter  VI. 

1.2  literature  Survey 

The  field  of  optimum  structural  design  has  been  exhaustively 
surveyed  frem  Galileo  to  the  present  in  review  papers  by  Wasiutynski 
and  Brandt  [2*0  and  Sheu  and  Prager  [22].  In  addition  a  review  by 
Gerard  [9],  which  is  particularly  applicable  to  aerospace  vehicles , 
contains  numerous  structural  references.  For  extensive  bibliographies 
covering  all  aspects  of  optimal  design  and  evaluations  of  the  current 
state  of  the  art,  the  reader  is  referred  to  these  articles. 

In  recent  years,  a  nvEber  of  investigators  have  applied  the 
methods  of  nonlinear  programming  to  optimum  structural  design.  Schmlt , 
Richer,  and  Morrow  [20]  solved  the  problem  of  integrally  stiffened 
waffle  plates  using  a  method  of  alternate  steps.  In  this  method 
steepest  descent  moves  are  made  until  a  constraint  is  violated.  The 


step  is  then  adjusted  until  the  design  is  on  or  near  a  constraint. 

This  point  is  considered  to  be  bound  and  an  alternate  step  to  a  free 
(unconstrained)  point  is  sought.  The  technique  uses  a  randan  number 
generator  to  create  proposed  alternate  step  designs  vhich  have  the 
same  weight  as  the  current  bouDd  design  and  which  do  not  violate  addi¬ 
tional  constraints.  Gellatly  and  Gallagher  [73  presented  the  theo¬ 
retical  basis  for  a  design  procedure  which  includes  steepest  descent 
and  alternate  step  moves  similar  to  those  of  Schait,  et.  al.  [20]. 

They  apply  this  procedure  to  the  design  of  plane  trusses  and  stiffened 
panels  [83.  The  objective  function  is  linear  in  both  cases,  since 
only  one  design  variable  per  structural  element  is  considered. 

Others  have  also  shown  that  optimun  structural  design  can  be 
formulated  as  a  nonlinear  programming  problem.  Pope  [153  solved  plane 
trusses  by  reducing  the  resulting  nonlinear  programming  problem  to  a 
sequence  of  linear  programming  problems,  and  Best  [2]  suggested  the 
use  of  a  gradient  projection  technique.  Recently  Brown  and  Ang  [33 
applied  the  gradient  projection  method  to  the  elastic  design  of  WF 
steel  frames.  A  problem  :lnvolving  four  design  variables  and  eleven 
constraints  was  solved  and  forward  finite  difference  approximations 
were  used  to  evaluate  the  derivatives  required. 

A  somewhat  different  approach  was  suggested  by  Scbmit  and  Fox 
[213.  Instead  of  developing  designs,  analysing  them,  and  using  infor¬ 
mation  from  the  objective  function  and  the  constraints  to  predict 
improved  designs,  a  special  function,  i|>,  is  constructed  and  minimised. 
The  \f>  function  contains  among  other  terms  the  weight  (or  other 
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quantity  to  be  minimized)  and  an  estimate  of  the  minimum  weight.  The 
function  is  no  structured  that  when  «  0  all  constraints  and  the 
equations  of  equilibria  and  compatibility  are  satisfied.  In  addition 
the  weight  is  decreased  to  the  estimated  minimum.  The  estimate  of 
the  minimum  weight  is  then  reduced  and  the  process  is  repeated  until 
the  weight  cannot  be  further  decreased.  At  each  step  the  solution 
i|  =  0  is  found  by  methods  for  unconstrained  minimization. 

A  method  of  nonlinear  programming  known  as  the  sequential  uncon¬ 
strained  minimization  technique  (SUMT),  has  also  received  some  atten¬ 
tion.  Hicholls  Ci43  used  the  method  to  solve  sane  plane  truss  prob¬ 
lems,  but  had  limited  success.  The  lack  of  success  was  attributed  to 
the  use  of  a  first  order  gradient  technique  for  the  unconstrained 
minimizations.  Brown  and  Ang  [4]  have  used  SUMT  to  obtain  starting 
values  for  the  gradient  projection  method  and  have  noted  the  possi¬ 
bility  of  solving  structural  design  problems  using  SUMT  instead  of 
other  programming  techniques.  A  recent  book  by  Bracken  and  McCormick 
[1]  discusses  examples  of  nonlinear  programming  problems  to  vhich  SUMT 
has  been  successfully  applied.  Among  these  is  the  design  of  a  verti¬ 
cally  corrugated  transverse  bulkhead  for  an  oil  tanker.  The  design 
is  specified  by  six  design  variables  and  must  satisfy  sixteen  con¬ 
straints.  The  objective  function  is  nonlinear  and  there  are  both 
linear  and  nonlinear  constraints.  Derivatives  of  the  constraint 
functions  were  available  without  approximation. 

The  programming  problems  associated  with  optimum  design  sire 
often  nonconvex.  In  general,  therefore,  a  local  minimum  is  obtained. 


Usually  widely  separated  starting  values  are  chosen  and  if  the  method 
converges  to  the  same  point  each  time,  this  point  is  taken  as  the 
global  minimum.  Toakley  [233  discusses  the  problem  of  the  global 
optimum  for  statically  determinate  plane  trusses.  A  change  of  vari¬ 
able  is  introduced  to  obtain  a  convex  nonlinear  programming  problem 
which  is  then  solved  by  reducing  it  to  a  sequence  of  linear  program¬ 
ming  problems . 

The  literature  concerning  mathematical  programming  has  also  been 
voluminous  in  the  past  several  years.  Recent  bookB  by  Wilde  and 
Beightler  [25],  and  Saaty  and  Bram  [193  discuss  a  wide  variety  of 
optimization  techniques.  The  former  treats  all  types  of  problems 
while  the  latter  has  an  extensive  section  on  the  methods  of  nonlinear 
programming.  A  book  by  Fiacco  and  McCormick  [53  presents  the  theo¬ 
retical  basis  of  the  sequential  unconstrained  minimization  technique. 
In  addition,  these  hooks  serve  as  a  review  and  bibliography  of  all  of 
the  major  contributions  to  the  field. 

Fletcher’s  review  [6]  of  unconstrained  minimization  techniques, 
which  do  not  require  derivatives,  is  particularly  pertinent  to  this 
research.  It  suggests  that,  of  the  methods  available,  the  one  due  to 
Powell  [16]  based  on  conjugate  directions  is  the  most  effective, 
Wortman  [253  has  written  a  program  which  combines  SUMT  and  Powell's 
method.  He  reported  extreme  sensitivity  to  the  starting  values  and 
proposed  that  the  technique  could  be  used  to  determine  binding  con¬ 
straints.  These  constraints  could  then  be  used  to  decrease  the 
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dimension  of  the  optimisation  problem  and  SUMT  could  then  be  reapplied. 
No  results  of  this  procedure  were  included  in  the  report. 

The  growing  literature  available  on  SUMT  [5  and  ll  indicates  that 
the  method  is  most  successful  when  derivatives  of  both  the  objective 
and  constraint  functions  can  be  computed  exactly.  One  of  the  objects 
of  this  dissertation  is  to  consider  the  applicability  of  SUMT  to  struc¬ 
tural  design  problems  when  derivatives  are  not  available. 

As  stated  previously,  almost  ail  design  problems  have  an  objec¬ 
tive  function,  design  variables,  state  variables,  a  set  of  state  equa¬ 
tions,  state  constraints,  and  design  constraints.  The  nonlinear  pro¬ 
gramming  approaches  cited  above  (except  for  Schmit  and  Fox  [21]),  use 
the  ntate  equations  indirectly  to  express  all  of  the  constraints  as 
functions  of  the  optimization  variables.  The  method  of  constrained 
steepest  descent  with  state  equations,  to  be  developed  in  Chapter  IV, 
differs  from  these  methods  in  that  the  state  equations  and  state  vari¬ 
able  constraints  appear  explicitly  in  the  mathematical  programming 
problem.  This  provides  a  very  natural  matching  of  the  essential  fea¬ 
tures  of  the  design  problem  and  the  method  used  to  obtain  its  solution. 
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CHAPTER  II 

FORMULATION  OF  THE  PROBLEM 

2.1  The  Problem 

This  chapter  presents  a  general  formulation  for  the  optimal 
design  of  minimus  weight  spatial  structures,  subjected  to  geometrical 
as  well  as  stress  constraints.  The  geometrical  constraints  are 
restrictions  placed  on  the  dimensions  of  the  structure,  whereas  the 
stress  constraints  correspond  to  a  failure  criterion  established  by 
one  of  the  failure  theories.  The  elements  of  such  a  structure  may 
undergo  extension,  bending,  twisting,  and  shear  deformation.  The 
formulation  of  the  design  problem  in  this  chapter  is  general,  but  for 
clarity  in  presentation,  some  examples  which  are  characteristic  of 
the  general  class  of  problems  will  be  used  in  the  development . 

Consider  a  structure  which  consists  of  rectangular,  hollow  beams. 
These  beams  are  mutually  perpendicular  to  each  other  and  are  Joined 
at  their  ends  by  rigid  Joints  so  that  forces,  bending  moments,  and 
twisting  moments  are  transmitted  from  member  to  member.  The  structure 
is  loaded  by  concentrated  forces  located  arbitrarily  along  each  of 
the  members.  The  frame  shown  in  Figure  2.1  is  an  example  of  this 
type  of  structure;  Figure  2.2  shows  a  typical  member  cross-section. 
Bending  and  twisting  effects  in  each  member  will  be  considered,  but 
shear  due  to  transverse  loads  and  axial  deformation  of  the  members 
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will  be  neglected.  The  thickness,  t,  of  the  bollow  sections  is 
constant  and  all  members  are  assumed  to  be  made  of  the  same  linearly 


elastic  material. 

The  design  variables  in  this  problem  are  the  thickness,  t, 
width,  b,  and  height,  h,  of  the  individual  members  (see  Figure  2.2). 
These  are  to  be  chosen  so  that  the  objective  function  (the  weight  of 
the  structure)  is  a  minimum  and  the  constraints  on  stress  and  geometry 
are  satisfied.  The  state  variables,  describing  the  behavior  of  the 
structure  for  a  given  design,  are  the  stresses  at  critical  points  in 
the  frame  and  the  deflections  at  the  joints.  The  state  equations 
which  determine  the  state  variables  are  the  equations  of  structural 
analysis  which  are  discussed  in  the  next  section. 

2.2  The  State  Equations 

The  state  variables  are  determined  from  a  set  of  matrix  equa¬ 
tions  which  are  derived  from  energy  principles.  Matrix  methods  for 
the  analysis  of  structures  have  received  much  attention  in  recent 
literature.  Books  devoted  to  this  subject  include  those  by  Ruben- 
svein  [lfi],  Przemieniecki  [17],  and  Zienkievicz  [273.  The  analysis 
of  the  structures  investigated  in  this  research  follows  a  method 
outlined  by  Langhaar  [11].  When  this  method  is  written  in  matrix 
notation,  the  resulting,  equations  are  of  the  same  form  as  those  used 
in  the  finite  element  and  matrix  displacement  methods  of  structural 
analysis.  Details  of  the  analysis  are  given  in  the  appendix;  the 
principal  features  of  the  method  are  as  follows. 


The  total  potential  energy  of  a  structure  and  the  »,ppli«>d  loads 


nay  he  written  as 


Y(< i)  ■  U(q)  <•  8(q) 


(2-1) 


where  U(q)  is  the  strain  energy  of  deformation,  Q(q)  is  the  potential 
energy  of  the  external  loads,  and  $  is  a  rector  whose  elasents  are  the 
cartesian  components  of  the  independent  kinematically  admissible 
angular  and  linear  joint  displacements.  Tor  sufficiently  small  dis¬ 
placements  the  internal  strain  energy  U{q)  is  a  positive  definite 
quadratic  form  in  the  generalised  coordinates  q.. 


inn 

U*2  1  2  aij*iV 

i«l  J«1  J  0 


(2-2) 


It  is  shown  in  Lenghsar's  hook  [11)  that  the  required  conditions  for 
equilibria*  are 


i  —  1,  2,  n. 


(2-3) 


where  u  is  the  number  of  generalised  coordinates.  This  principle  of 


stationary  potential  energy  applied  to  equation  (2-1)  yields 
|L.+  |S-»°  i-1,  2,  ...,  n. 


(2-4) 


Using  the  definition  of  the  components  of  generalised  force  given 


(2-5) 


the  following  set  o.  . qvations  is  obtained 


A  *****  m  Pi? 


i  »  1,  2, 


In  matrix  notation  these  equations  nay  he  expressed  as 
AQ  *  ? 


(2-6) 


(2-7) 


*  jSUKmtf 
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which  corresponds  directly  to  the  basic  equation  of  the  matrix  dis¬ 
placement  method  of  structural  analysis  [171. 

The  elements  of  matrix  A  are  known  functions  of  the  design  vari¬ 
ables  and  the  elements  of  vector  P  are  functions  of  the  external  loads. 
Equation  (2-7)  may  be  solved  for  the  joint  displacements  Q.  In  the 
appendix  it  is  shown  that  the  end  reactions  and  the  member  stresses 
may  also  be  written  as  matrix  equations.  For  end  reactions  the  fol¬ 
lowing  equation  is  obtained 

M  =  BQ  +  F,  (2-8) 

in  which  M  is  the  vector  of  end  reactions ,  B  is  a  matrix  whose  ele¬ 
ments  are  functions  of  the  design  variables,  Q  is  the  vector  of  Joint 
displacements,  and  P  is  e  vector  of  the  "fixed-end"  mounts  and  loads. 
The  stresses  may  be  computed  from  the  following  equation 

S  =  CM  +  ?  (2-9) 

where  S  is  the  vector  of  stresses  calculated  at  critical  points  in 
the  structure  (see  Appendix),  C  is  a  matrix  whose  elements  are  func¬ 
tions  of  the  design  variables,  M  is  the  vector  of  end  reactions  *  and 
P  is  a  vector  whose  elements  are  functions  both  of  the  applied  loads 
and  the  design  variables.  By  combining  equations  (2-8)  and  (2-9), 
the  following  equation  is  obtained 

S  a  CBQ  -v  CF  +  P.  (2-10) 

The  analysis  of  the  structure  may  then  be  summarized  by  the  fol¬ 
lowing  equations 

AQ  =  P 
M  =  BQ  +  F 

and  S  =  CBQ  +  CF  +  P.  (2-11) 


\ 


i 
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The  fundamental  equations  of  the  finite  element  technique  are  given 
by  Zienkievicz  [273  and  are  of  the  same  form  as  equations  (2-ll). 
Conceptually  the  optimization  procedures  investigated  in  this 
research  can  then  be  applied  to  the  broad  classes  of  problems  which 
can  be  analyzed  by  matrix  displacement  and  finite  element  techniques. 
Equations  (2-11)  apply  to  linearly  elastic  structures  under  the 
action  of  conservative  external  loads. 

Physically  the  behavior  of  the  frame  is  most  often  described  by 
the  joint  deflections  and  stresses.  Therefore  it  is  natural  to 
define  these  variables  as  the  state  variables,  Z, 


Z  * 


(2-12) 


t  J 

Given  the  external  loads  and  a  set  of  design  variables,  the  first 
and  third  equations  of  set  (2-11)  can  then  be  combined  to  yield  one 
matrix  equation  which  determines  Z.  This  equation  is  given  by 


AZ  -  P 


* 


(2-13) 


where 


A  = 


and 


J 


[CF+PJ  . 

The  analysis  of  the  structure  (determination  of  the  state  variables) 
*  &  therefore  reduced  to  the  solution  of  a  single  set  of  simultaneous 


equations . 
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2.3  The  Objective  Function 

Since  all  members  are  assumed  to  be  made  of  the  3araa  material, 
the  weight  of  the  frame  is  the  product  of  a  single  weight  density, 
and  the  volume  of  the  frame.  Therefore  the  problem  of  determining 
the  minimum  weight  is  equivalent  to  that  of  determining  the  minimum 
volume.  The  volume  of  the  ith  member  of  the  frame  is  given  by 

vi  ■  tVi  -  <bi  -  2V(hi  -  2tinci 


or 


V.  =  [2tlbi  ♦  Sit A  -  UtfKj 


(2-14 ) 


1U 


where  t^ ,  b^ ,  and  h^  are  defined  in  Figure  2.2,  and  (b  is  the  length 
of  the  member.  The  design  variables  are  systematically  assigned  as 
follows ,  m  being  in  this  case  the  number  of  frame  members , 

xi  =  ti 

X2  “bl 

*3  *  hl 


x3i-a  1  h 


x31-l  =  bi 


x3i  ‘  hi 


i  —  1,  . . . ,  m« 


Using  this  notation,  equation  (2-lU)  can  be  written  in  terms  of  tbe 
design  variables: 

Vi(x)  =  C2x3i_2X3i_i  +  2x3i-2X3i  “  4x‘"3i-2-1Ci*  ^2“15^ 
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The  volume  of  the  frame  is  obtained  by  summing  the  volumes  of  the 
individual  members  and  is  a  quadratic  in  the  design  variables . 


m 

V(x)  =  Z  V  (x) 
i=l 


(2-1 6) 


2.1*  The  Constraints 

Structural  failure  is  said  to  occur  when  the  state  variables  or 
functions  of  the  state  variables  exceed  certain  preassigned  limits. 

In  space  frames  failure  may  be  due  to  either  excessive  stresses  or 
deflections.  In  the  typical  structures  being  considered,  failure  is 
assumed  to  occur  when  the  combination  of  normal,  stress  due  tc  bending 
and  shear  stress  due  to  torsion  exceeds  a  maximum  limit.  The  stresses 
were  combined  according  to  the  distortion  energy  theory  of  failure. 

The  failure  stress  a4  is  confuted  as  follows: 


o2  +  Sr1 


(2-17) 


where  a  is  the  normal  stress  and  t  is  the  shearing  stress  at  a  criti¬ 
cal  point  of  the  structure.  The  statement  that  the  failure  stress 
must  not  exceed  a  specified  maximum  is  written 

/a2  -  3tz  -  a  <0  (2-l8) 

max 

and  the  general  form  of  a  state  variable  constraint  is 

iKz)±0.  (2-19) 

The  choice  of  design  variables  may  also  be  limited.  These  limits 
may  occur  because  of  space  restrictions  or  result  from  limits  imposed 


by  current  manufacturing  facilities.  Constraints  of  this  type  are 
introduced  into  the  structural  design  problem  of  interest  by  requiring 


that  the  design  variables  satisfy  the  following  inequality: 


X  .  <  X  <  X 

min  —  —  max 


(2-20) 


This  inequality  is  equivalent  to  two  constraint  relationships,  namely. 


x  .  -  x  <  0 

min  — 


(2-21) 


x  -  x  <0 
max  — 


(2-22) 


The  general  form  of  these  constraints  is 
4>(x)  <  0. 


(2-23) 


2.5  The  Mathematical  Programming  Problem 

The  state  equations,  the  objective  function,  and  the  constraints 
associated  with  the  optimal  structural  design  problem  have  been 
defined  in  the  previous  sections  of  this  chapter.  These  functions 
may  be  used  to  define  the  following  mathematical  programming  problem: 


minimize  Y(x) 


(2-2U) 


suoject  to 


^(x)  <_  0, 
^(z)  <  0, 
h^Zjx)  =  0, 


i  •'  1,  ...,  I 
i  *  1,  . . . ,  m 
i  ®  1,  k . 


The  equality  constraints  in  problem  (2-24)  are  the  state  equations 
(2-13).  The  mo3t  common  approach  to  this  problem  has  been  to  use 
the  state  equations  to  write  the  state  constraints  as  functions  of 
the  design  variables.  This  results  in  the  following  nonlinear 


programming  problem,  which  may  be  attacked  by  several  available  tech¬ 
niques  (see  Section  1.2), 

minimize  V(x)  (2-25) 

subject  to 

gj(x)  <  0  j  “  1,  ...»  m +1 

where 


A  design  which  satisfies  ill  of  the  impose  *  constraint  conditions  is 
called  a  feasible  design.  The  constraints  define  a  feasible  region, 
R,  in  design  space, 

R  =  (x:  g(x)  £  0}.  (2-26) 

If  the  nonlinear  programming  problem  (2-25)  is  convex,  then  it  can  be 
shown  that  the  solution  attained  is  a  global  optimum.  It  will  now  be 
shewn  that  the  structural  optimization  problem  formulated  in  this 
chapter  is  nonconvex. 

The  convexity  of  sets  and  functions  may  be  investigated  by  using 
three  theorems  which  are  proven  in  (.19]  . 

Theorem  2.1  The  set  of  points  R  which  satisfy  a 
constraint  g(x)  <_  0,  where  g(x)  is  n  convex  func¬ 
tion  is  a  convex  set. 

Theorem  2.2  The  intersection  R  of  a  family  F  of 
convex  sets  is  a  convex  set. 


wmm 
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Theorem  2.3  If  V(x)  is  a  twice  differentiable 

function  in  an  open  convex  set  R,  it  is  convex 

in  R  if  and  only  if  the  quadratic  fora 
n  n  a2v 

i=l  j=l  8xidxj  1  J 

is  positive-semidefinite  for  every  point  x  in  R. 


Because  the  state  variable  constraints  are  complicated  functions 
of  the  design  variables,  it  is  difficult  to  apply  Theorem  2.3  to 
these  functions.  Therefore  one  cannot  immediately  determine  whether 
or  not  the  feasible  region,  R»  i.J  a  convex  set.  Consequently,  the 
convexity  of  the  objective  function  will  be  investigated  by  applying 
the  criterion  of  Theorem  2.3  assuming  that  the  region  R  satisfies  the 
conditions  of  the  theorem. 

Consider  the  volume  of  the  three  member  frame  shown  in  Figure  2.1 

V(x)  =  (2x^2  *  2xix3  ~  ^)Ci 
+  tex^x*.  +  2xux6  -  Uxjj)C2 
*  (2x^Xg  +  2x^xQ  “  *tx^)C^.  (2-27) 

The  matrix  of  second  derivatives  required  in  Theorem  2.3  is 


The  analysis  of  quadratic  forms,  outlined  in  Langhaar's  book  [11], 
may  be  applied  to  show  that  the  quadratic  form  based  on  this  matrix 
is  indefinite.  Consequently,  by  Theorem  2.3,  V(x)  is  a  nonconvex 
function  when  the  region  R  is  convex.  Therefore  problem  (2-25)  is 
a  nonconvex  programming  problem. 

The  most  that  can  be  guaranteed  for  a  nonconvex  problem  is  that 
its  solution  is  a  local  minimum.  The  most  often  used  technique  for 
seeking  the  global  optimum  in  this  situation  is  to  attempt  to  deter¬ 
mine  all  of  the  local  minima  by  starting  from  many  widely  separated 
initial  points. 

In  summary  the  optimal  structural  design  problem  has  been  formu¬ 
lated  as  a  mathematical  programming  problem  in  two  different  ways. 

In  one  case  the  constraints  are  considered  to  be  functions  of  the 
design  variables  alone.  The  sequential  unconstrained  minimization 
technique  will  be  applied  to  this  problem  in  Chapter  III.  In  Chapter 
IV  the  constrained  steepest  descent  method  with  state  equations,  will 


SSST'^PSPI 


■f^rsrr.-v&Sri  ssr  >, 


be  applied  to  problem  (2-24) .  This  formulation  considers  constraints 
which  are  functions  of  both  the  state  variables  and  the  design  vari¬ 
ables,  and  uses  the  state  equations  directly  in  the  solution. 
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CHAPTER  III 

THE  SEQUENTIAL  UNCONSTRAINED  MINIMIZATION  TECHNIQUE 


gi(x)  <  0  i  *  1,  2,  ....  m+l 

The  constraints  are  complicated  nonlinear  functions  of  the  design 
variables.  Because  of  this,  the  direct  handling  of  these  con¬ 
straints  in  the  nonlinear  programing  problem  can  be  difficult. 

The  sequential  unconstrained  minimization  technique  (SUMT)  handles 
these  constraints  indirectly  and  has  been  used  successfully  on 

problems  of  this  type.  The  theoretical  basis  of  this  technique, 

as  veil  as  helpful  suggestions  for  computation,  may  be  found  in 

the  recent  book  by  Fiacco  and  McCormick  [53.  There  are  several 

versions  of  SUMT  which  may  be  applied  to  given  problems  depending 

upon  the  nature  of  the  objective  function  and  the  constraints.  For 

problem  the  interior  method  should  be  used.  In  this  method 

the  quest  for  a  minimum  is  always  carried  out  within  the  feasible 

region  and  strict  equality  constraints  arc  not  allowed. 

The  method  is  applied  by  augmenting  the  objective  function  to 

define  a  new  function  U(x,rk), 
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U(x,rk)  .  V(x)  +  S(rk)l(x).  (3-2) 

The  number  rk  is  always  positive.  l(x)  is  a  scalar  function  of  the 
design  variables  and  is  continuous  in  the  feasible  region  R,  vhere  R 
is  defined  by 

R  *  (x:  g^(x)  <  0,  i  ■  1,  2,  ...,  m+£). 

I(x)  also  has  the  property  that  if  (x^}  is  any  infinite  scqueaee  of 
points  in  R  converging  to  Xg  such  that  g^Xg)  *  0  for  at  least  one 
it  then 

lim  Kx.  )  *  +  •. 

S(rk)  is  a  scalar-valued  function  of  the  single  variable  r^,  and  has 

the  following  properties.  If  r1  >  rg  >  0,  then  S(r1>  >  S(rg)  >  0, 

and  if  {rfe}  is  an  infinite  sequence  of  points  such  that 

lim  rk  *  0 
k-*» 

then 

,  lim  S(rt)  *  0. 

These  properties  are  basic  to  the  convergence  proofs  which  may  be 
found  in  Fiacco  and  McCormick  [53. 

The  most  com non  forms  for  S(?k)  sad  I(x)  are  the  fallowing: 


S(rk)  *  r 

(3-3) 

m*Z  , 

I(x)  *  -  Z 
i»l 

'  (3-fc) 

Using  these,  the  function  U(x,rk>  becomes 
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Ml  ^ 

U(x,r.  )  ■  V(x)  -  r  I  r-T-y  *  *  °*  ^3“^) 

K  i*l  ®iv/  A 

The  penal fry  term  pj  o  video  &  steep  gradient  at  the  boundary  ov  the 
feasible  region  no  that  when  the  function  U(x„rk)  is  minimised  the 
search  is  kept  in  the  feasible  region  R.  When  g^x)  _>  0,  U'x.r^)  ■  K, 
where  K  is  a  large  positive  constant.  This  avoids  difficulties  which 
may  arise  in  the  unconstrained  minimisation  if  a  constraint  is  inad¬ 
vertently  violated. 

The  algorithm  for  finding  the  constrained  minimum  of  a  function 
using  interior  SUMT  is  as  follows. 

(1)  Choose  rx  and  an  estimate  of  the  minimum  x^  interior  to  R. 

(2)  Determine  the  unconstrained  minimum  x(r^)  of  U(x,:r]L)  in  the 
feasible  region  R. 

(3)  Use  x(r1)  as  a  new  starting  point  to  determine  the  minimum 
x(r2)  of  UfXfig)  '•'here  r±  >  rg  >  0. 

(k)  Continue  the  process,  finding  the  local  rsiniunrr  of  U(x,r^) 
starting  from  x(rfc  ^  for  a  strictly  aonotonically  decreasing 

sequence  (r^). 

Fiacco  and  McCormick  show  that  if  V(x)  and  g<(x),  i  ■  1,  2,  ...,  m*£ 
are  continuous,  if  l{x)  and  S(r k>  satisfy  the  conditions  stated  previ¬ 
ously,  if  there  exists  a  relative  minimum  x  in  R  such  that  f(x>  >  f(x) 
for  all  x  ^  x  in  some  neighborhood  of  x,  and  x  is  not  an  isolated 
point  in  R,  and  if  trk>  i«  *  sonotone  decreasing  sequence  which  con¬ 
verges  to  aero,  then  the  sequence  of  local  unconstrained  minima 


I  . 


K-  .  ( 


2k 


obtained  in  the  above  procedure  exists  and  converges  to  a  local 
minimum  of  problem  (3-l). 

The  use  of  SUifF  is  appealing  in  that  existing  unconstrained 
minimi nation  techniques  can  be  used  to  obtain  the  minimum  for  each 
U(x,rk).  The  method,  therefore,  provides  a  means  of  reducing  the 
value  of  the  objective  function  without  violating  tho  constraints. 
Furthermore  it  is  not  necessary  to  involve  the  constraints  directly 
or  to  move  along  them  toward  the  constrained  minimum  of  the  objective 
function.  This  is  a  distinct  advantage  when  the  constraints  are  com¬ 
plicated  nonlinear  functions  of  the  design  variables  as  they  are  in 
the  class  of  problems  under  consideration. 

In  the  next  section  an  unconstrained  minimisation  technique 
which  does  not  require  derivatives  will  be  discussed.  Acceleration 
of  convergence  by  extrapolation  will  be  discussed  in  Section  3.3, 
and  other  computational  considerations  will  be  presented  in  Section 
3.4,  An  evaluation  of  the  method  as  applied  to  the  optimum  struc¬ 
tural  design  problem  at:,  be  found  in  Chapter  V. 


3.2  Finding  an  Unconstrained  Minim  via 

As  stated  in  Chapter  I,  one  of  the  goals  of  this  research  is  to 
investigate  techniques  that  will  he  effective  when  derivatives  of  the 
objective  and  constraint  functions  are-  cumbersome  or  impossible  to 
calculate.  The  method  used  to  determine  the  unconstrained  minimus  of 
U(x,r,  )  in  SUMT  must  therefore  be  able  to  find  a  mini  mum  using  rally 


function  values. 


The  earlier  techniques  such  as  tabulation,  random  search,  or 
improving  one  variable  at  a  time  are  basically  inefficient  and  unre¬ 
liable.  The  first  improvements  of  these  methods  were  based  on  ad  hoc 
procedures.  Only  recently  efficient  techniques  have  been  introduced 
which  are  based  on  successive  minimizations  along  conjugate  directions. 

In  his  review  of  minimization  techniques  that  do  not  require  the 
calculation  of  derivatives,  Fletcher  [6]  suggests  that  of  the  methods 
available,  Foweli ' s  [16]  is  the  most  satisfactory.  Powell's  method 
requires  fever  function  evaluations  than  other  techniques  and  has  tne 
advantage  of  quadratic  convergence  near  the  minimum.  The  method  will 
find  the  minimum  of  a  quadratic  in  a  finite  number  of  steps  and  con¬ 
verges  to  the  minimum  from  an  unfavorable  starting  point  more  effi¬ 
ciently  than  other  available  algorithms. 

Powell's  method  is  based  on  the  minimization  of  a  quadratic, 

f(x)  =  x^[A]x  +  bx  +  c.  (3-6) 

Directions  and  £,>  are  said  no  he  conjugate  with  respect  to  A  if 

^[AK2  =  0.  (3-7) 

Each  iteration  starts  from  the  best  previous  estimate  to  the  minimum 
Pq.  Successive  searcnes  for  a  minimum  are  made  along  a  set  of 
linearly  independent  directions,  £^,  •••»  5n*  These  directions 

are  initially  chosen  as  the  coordinate  directions  ro  that  the  first 
iteration  is  identical  to  that  of  changing  one  parameter  at  a  time. 
Each  iteration  generates  a  new  direction  S,  and  for  the  second  iter¬ 
ation  the  set  of  linearly  independent  directions  is  chosen  to  be 
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r  r  The  new  direction  £  is  chosen  so  that  if  a  quad- 

2*  3*  n 

ratic  is  being  minimized  all  the  directions  vill  he  conjugate  after  n 
iterations.  Povell  proves  that,  as  a  direct  consequence,  the  exact 
minimum  of  the  quadratic  is  found.  One  iteration  of  the  basic  method 
is  as  follows. 

(i)  For  r  *  l,2,...,n  calculate  Xy  so  that  f(pr-j_  + 
is  a  minimum  and  define  pf  =  Pr_i  +  ^r^r* 

(ii)  For  r  =  1,2,..., n-1  replace  £r  by  5^. 

(iii)  Replace  ?n  by  (pn  -  pQ). 

(iv)  Choose  X  so  that  f(pQ  +  X(pn  -  PQ))  is  a  minimum 
and  replace  pQ  by  pn  +  X{pQ  -  Pq}- 

The  basic  procedure  may  be  unstable  for  non-quadratic  functions 
because  it  tends  to  choose  nearly  dependent  directions.  Powell  shows 
that  this  difficulty  can  be  eliminated  by  using  the  following  proce¬ 
dure  which  allows  a  direction  other  than  ^  to  be  discarded  and  under 
some  conditions  uses  the  old  set  of  linearly  independent  directions 
again. 

(i)  For  r  =  l,2,...,n  calculate  so  that  f(pr_^  +  ^r?r) 
is  a  minimum  and  define  pr  =  Pr_^  +  *r^r’ 

(ii)  Find  the  integer  m,  1  £  m  <_  n,  so  that  ff(pm_^)  ~ 
is  a  maximum,  and  define  A  =  f(pjn_1)  "  T(pm)» 

(iii)  Calculate  f 3  =  f(2pn  -  pQ)  and  define  =  f(pQ)  and 

t2  -  f(p„>- 

(iv)  If  either  f 3  >.  t±  and/or 

(fx  -  2fg  +  f3)*(f1  -  *2  ~  A)2  >  |  A  (^  -  f3)2 
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use  the  old  directions  5^,  5^,  ...»  Cn  for  the  next  iteration 
and  use  pfl  for  the  next  pQ,  otherwise 

(v)  defining  £  =  (p  -  prt),  calculate  X  so  that  f(p  +  X£)  is  a 

n  u  n 

minimum,  use  ^  «2,  ....  ^  ^+1>  Sm+2,  ....  C  as 

the  directions  and  pfi  +  X£  as  the  starting  point  for  the  next 
iteration . 

Complete  theoretical  justification  of  the  above  algorithm  is  given  in 
Powell's  paper.  One  effect  of  step  (v)  above  is  that  one  of  the  previ¬ 
ously  determined  conjugate  directions  may  be  thrown  away.  In  this 
case  the  minimum  of  a  quadratic  will  require  more  than  n  iterations. 
However,  Powell  shows  that  this  procedure  ensures  that  the  rate  of 
convergence  will  always  be  reasonable,  therefore  making  the  modifica¬ 
tion  valuable.  In  fact,  Powell  found  it  to  be  essential  when  mini¬ 
mizing  a  function  of  twenty  variables,  and  highly  desirable  for  func¬ 
tions  of  five  variables  or  more, 

Powell  suggests  a  very  safe  but  lengthy  convergence  criterion. 

This  procedure  was  not  used.  Instead,  in  order  to  decrease  the  num¬ 
ber  of  function  evaluations  required,  the  iterations  were  terminated 
when  the  results  from  two  successive  iterations  agreed  to  within  a 
specified  value,  e  .  According  to  Fletcher  [6],  only  the  most  diffi¬ 
cult  functions  require  the  more  stringent  convergence  criterion. 

If  Powell's  method  is  to  be  effective,  it  is  essential  that  an 
efficient  method  of  finding  a  minimum  along  a  given  direction  he 
available.  The  objective  function  surface  U(x,r^)  is  peculiar  in 
that  at  the  boundary  of  the  feasible  region  the  function  takes  on  a 


large  positive  value.  Taking  this  into  account  the  following  proce¬ 
dure  was  developed. 

(1)  Starting  from  p^  ^  use  a  Fibonacci  search  along 
the  direction  to  bracket  the  minimum  within 
the  feasible  region. 

(2)  Pass  a  quadratic  through  the  three  points  brack¬ 
eting  the  minimum  and  approximate  the  minimum  of 
the  function  along  the  line  by  the  minimum  of  the 
quadratic . 

(3)  Let  the  quadratic  be  determined  by  function  values 
corresponding  to  design  points  x^-  x2,  and  x^, 
and  let  the  minimum  of  the  quadratic  be  .  If 

l\in  ‘  xil  ‘  V  1  ‘  2’  3 

terminate.  Otherwise  retain  the  three  points 
which  bracket  the  minimum  and  reduce  the  interval 
of  uncertainty,  and  repeat  step  (2). 

The  above  algorithm  was  used  in  conjunction  with  Powell's  method  and 
provided  convergence  comparable  to  *v  reported  by  Powell  in  his 
paper. 

3.3  Extrapolation 

Fiaceo  and  McCormick  show  that  the  convergence  of  SUMT  can  be 
accelerated  by  extrapolation.  If  p  minima  of  U(x,r^)  have  been 
determined,  these  may  be  used  to  estimate  the  optimum  (the  minimum  of 


V(x) )  and  the  (p+l)st  minimum  of  U(x,rk). 


estimates  of  the  optimum  are  shown  to  converge  more  quickly 
to  the  solution  than  the  sequence  of  unconstrained  minima.  In  addi¬ 
tion  the  estimate  of  the  (p+l)st  minimum  can  "be  used  as  a  starting 
value  for  the  (p+-l)st  iteration  which  significantly  reduces  the 
amount  of  computation  required  to  find  the  optimum. 

After  p  minima  have  been  obtained,  estimates  of  order  p-1  can  be 
made.  In  practice,  estimates  are  seldom  made  beyond  order  three  due 


to  computer  storage  requirements  and  accuracy  considerations  such  as 
round-off  error.  The  experience  of  Fiacco  and  McCormick  is  that  even 


first  and  second  order  approximations  of  the  next  U(x,rk)  minimum  and 
the  optimum  significantly  accelerate  the  convergence. 

The  extrapolation  is  based  on  the  fact  theit  the  p  minima  which 


have  been  found,  x^,  Xg,  ...»  x^  which  correspond  to  r^,  r,,,  ...» 


may  he  expanded  in  terms  of  the  r^  as  follows.' 


P-1  . 

*k  =  1  ai*rk)J 

j=0  J 


k  =  1,  2,  . . . ,  p 


(3-8) 


where  are  n  component  vectors.  A  set  of  recursion  relations 
based  on  this  expansion  leads  to  the  following  equations  for  first 
and  second  order  estimates  of  the  optimum  and  the  (p+l)st  minimum  of 


the  U  function  when  r. 


r^/c,  (c  >  1). 


Hie  first  order  estimate  of  the  optimum  is, 


CX  -  X 


09) 


and  the  second  order  estimate  is  given  by 
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x*  = 


c3x  -  c2x  ,  -  ex  -  +  cx  _ 
_ 2 _ P-1  .  P-1  .  P-2 

(c2  -  l)(c  -  1) 


(3-10) 


The  first  order  estimate  of  the  (p+l)st  minimum  is 

x  +  cx  -  x  , 

x  «  _£ _ 2 _ £=1 

P+1 


(3-11) 


and  the  second  order  estimate  may  be  written 


c3x  +  c2x  -  c2x  ,  -  x 


P+1 


P-1  P-2 


(3-12) 


These  equations  were  used  to  accelerate  convergence  when  SUMT  was 
applied  to  the  structural  optimization  problem  being  considered.  The 
first  order  estimates  were  applied  after  two  iterations  had  been  com¬ 
pleted  and  the  second  order  estimates  were  used  thereafter.  The 
estimate  of  the  (p+l)st  minimum  was  used  as  a  starting  value  for  the 
(p+l)st  iteration  only  if  the  estimate  was  a  feasible  point. 


3.4  Parameter  Selection  and  Convergence  Criteria 

The  experience  of  Fiacco  and  McCormick  indicates  that  the  con¬ 
vergence  of  SUMT  is  not  greatly  effected  by  either  the  choice  of  r^  or 
by  the  choice  of  the  factor  C  by  which  r^  is  reduced  at  each  iteration. 
This  was  a-5-  '  confirmed  by  the  author’s  experience  when  SIMT  was 
applied  to  the  structured,  optimization  problem  where  values  r^  *»  1, 
and  C  «  4  were  chosen. 

When  the  differences  between  components  of  two  successive 

of  the  U(x,r.  )  function  were  less  than  a  preassigned  value,  c  ,  the 
X  sx 

process  was  terminated.  The  choice  of  this  preassigned  value  is  not 
arbitrary.  Jn  particular  it  cannot  be  less  than  the  value  used  as  a 


convergence  criterion  in  the  unconstrained  minimization.  Similarly, 
the  convergence  criterion  for  the  one  dimensional  minimization  must 
he  most  stringent  of  all. 


f 
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CHAPTER  IV 

CONSTRAINED  STEEPEST  DESCENT  WITH  STATE  EQUATIONS 

H.l  Introduction 

The  nonlinear  programming  technique  of  Chapter  III  uses  the 
state  equations  indirectly  to  determine  the  values  of  the  constraints 
at  a  given  point  in  design  space.  One  way  of  using  the  state  equa¬ 
tions  more  explicitly  in  the  nonlinear  programming  problem  in  to 
introduce  them  as  additional  constraints  on  the  solution,  and  then 
to  use  the  design  variables  plus  the  state  variables  as  the  inde¬ 
pendent  variables  in  the  problem.  Oftentimes,  however,  this  may  not 
lead  to  a  satisfactory  solution,  since  in  many  problems  the  number 
of  state  variables  is  large  compared  to  the  number  of  design  vari¬ 
ables.  The  dimension  of  the  resulting  programming  problem  become* 
very  large.  For  instance,  one  of  the  simplest  examples  discussec 
in  Chapter  V  has  six  design  variables  and  ten  state  variables,  which 
would  then  lead  to  a  nonlinear  programming  problem  of  dimension  six¬ 
teen.  If  the  derivatives  of  the  objective  function  and  the  con¬ 
straints  are  cumbersome  to  compute  or  are  otherwise  not  available, 
the  situation  is  even  more  acute  in  that  the  most  successful  minimi¬ 
zation  techniques,  which  do  not  require  derivatives,  seem  to  be 
limited  to  problems  of  about  twenty  variables  [6  &  16]. 


In  thds  chapter  a  method  is  introduced  vhich  exploits  the  state 
equations  by  using  them  directly  in  the  design  process.  The  general 
scheme  is  to  linearize  the  problem  in  the  neighborhood  of  a  candi¬ 
date  design  in  terms  of  the  design  and  state  variables,  and  then  to 
eliminate  direct  dependence  upon  the  state  variables  by  introducing 
the  linearized  state  equations.  This  results  in  a  nev  nonlinear  pro¬ 
gramming  problem  to  vhich  the  Kuhn-Tucker  necessary  conditions  are 
applied.  As  a  consequence  a  step  is  chosen  in  design  space  vhich 
reduces  the  objective  function  consistent  vith  the  constraints  and 
simultaneously  directs  tne  search  for  an  optimum  back  tovards  the 
feasible  region  if  any  of  the  constraints  have  been  violated.  The 
method  was  developed  by  Haug  [10]  who  generalized  ideas  introduced  by 
Mel'ts  [13].  Since  the  method  is  not  yet  available  in  the  literature, 
it  will  be  discussed  in  detail  in  the  following  section.  Section  4.3 
will  present  the  basic  algorithm  and  discuss  some  of  the  computational 
aspects  of  the  method.  Convergence  will  be  considered  in  Section  4.4 
and  an  expanded  algorithm  is  presented  in  Section  4.5. 

4.2  Description  of  the  Method 

In  Chapter  II  it  was  shown  that  when  the  state  equations  are 
included  in  the  formulation,  the  optimal  structural  design  problem 
may  be  stated  as  the  following  mathematical  programming  problem: 


ass  ssi* 
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minimize  V(x) 

subject  to 

^(x)  £  0  i  ■  1,  2,  ...»  m  (!*-l) 

hi(z,x)  a  0  i  *  1,  2,  ...»  k 

^(z)  £  0  i  =  1,  2,  ...»  t. 

With  minimal  additional  effort,  the  method  of  constrained  steepest 
descent  with  state  equations  (CSDS)  can  be  derived  for  a  more  general 
case  which  allows  the  objective  function  and  the  <j>  constraints  to 
depend  upon  the  state  variables,  z,  as  well  as  the  design  variables, 
x.  The  more  general  problem  may  be  stated  as  follows: 
minimize  f (z,x) 

subject  to 

^(z,x)  £0  i  =  1,  2,  ...»  m  (U-2) 

h^Zjx)  =0  i  *  1,  2,  ...,  k 

♦jU)  £  o  i  *  1,  2,  ...»  1. 

If  desired,  the  <|i  constraints  could  be  considered  as  a  subset  of  the 
4>  constraints.  If  the  mathematical  form  of  the  two  types  of  con¬ 
straints  is  sufficiently  similar,  this  will  provide  some  simplifica¬ 
tion  in  the  computational  algorithm.  If  the  forms  of  the  constraints 
are  considerably  different,  sixne  computational  advantage  may  be 
gained  by  considering  them  separately.  For  example,  if  the  $  con¬ 
straints  are  very  simple  and  the  constraints  complicated,  then  it 
may  be  advantageous  to  compute  the  derivatives  of  the  <p  constraints 
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directly  but  to  approximate  the  derivatives  of  the  i>  constraints. 

In  this  section  the  derivation  for  the  design  problem  (**-2)  is  pre¬ 
sented 

ThxS  chapter  makes  use  of  matrix  calculus  notation,  A  function 
g  of  the  vector  x  will  be  defined  as  follows: 

gx(x) 


g(x) 


where 


x  - 


*1 


(x) 


The  derivative  of  g(x)  with  respect  to  x  is  a  matrix , 


3x 


!Sl 

3x 


Jj 


of  dimension  mxn.  A  small  change  in  x  will  be  denoted  by 

6xl 


<5x 


<5x- 


6x 


3T 


H  =  f(  z,x)  +  X^hCz.x)  +  +  vTiJ>, 

then  3H/3x  =»  0 

3H/3z  »  0 

*  0  if  <Ci(z,x)  <  0 
and  vi  «  0  if  ^(z)  <  0 


(4-3) 

(4-4) 

(4-5) 


at  x  “  x*  and  z  *  z*. 

One  could  attempt  to  use  equations  (4-4)  and  ( U— 5 )  of  the  theorem  to 
construct  solutions  of  problem  ( 4-2 ) ,  but  it  is  very  difficult  to 
determine  which  of  the  inequalities  are  strict  equalities  at  the  solu¬ 
tion  when  the  number  of  constraints  is  large.  Instead,  the  theorem 
will  be  used  in  this  section  to  develop  a  direct  method  of  solving 
the  optimal  design  problem. 

Let  x  be  an  estimate  of  the  solution  to  problem  (4-2)  and  z  be 
the  state  variables  associated  vith  this  solution.  The  purpose  of 
the  method  is  to  determine  a  small  change  in  x,  ox,  which  will 
decrease  the  objective  function  and  satisfy  the  constraints.  To 
obtain  this  goal,  the  problem  is  first  line  crimed  about  the  estimate 
to  the  solution  (x,z)  where  h(z,x)  *  0.  When  this  is  done,  the  fol¬ 
lowing  first  order  changes  in  the  functions  of  problems  (4-2)  result. 


6f  ■ 

|X6z 

3z 

,  + 

3f  . 
37  6x 

(4-6) 

3*. 

_  l 

3z 

5z 

+  sr 

Sx 

(4-7) 

3h. 

3h. 

6hi 

_  i 
~  3z 

oz 

♦sr 

6x 

(4-8) 
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3z 


6  z 


(M) 


It  should  he  remembered  that  the  partial  derivatives  in  equations 
(4-6)  through  (4-9)  are  evaluated  at  the  estimate  (x,z)  of  the  solu¬ 
tion.  The  change  in  the  design  variables  and  the  resulting  change  in 
the  state  variables  must  satisfy  h(z  +  6z,x  +  6x)  «  0.  therefore 
it  is  required  that  $h  *  0.  If  at  the  estimate  of  the  solution  the 
inequality  constraints  are  violated,  then  changes  6<J>^  *  A+^,  and 
<5^  *  Ai^i  are  requested  such  that  these  constraints  will  be  satisfied 
•within  the  linear  formulation  of  the  problem.  If  $^(z,x)  »  0  and 
^(z)  *  0  then  it  is  required  that  6+,  *  0  and  6^  ■  0.  In  addition 
the  accuracy  of  the  linear  approximation  must  be  guaranteed  by  insur¬ 
ing  that  the  step  size  remains  small.  Therefore  it  is  required  that 


6x  wdx  «  £2 


(4-10) 


where  5  is  a  small  number  and  w  is  a  positive  definite  weighting 
matrix.  The  linearized  version  of  the  problem  may  then  be  written  as: 

minimize  ~  6z  +  ~  Sx  (4-11) 

oZ  oX 


subject  to 


||  6z  +  ||-  6x  =  A* 


(4-12) 


3h  ,  .  ah  . 
(I*1’!!*1"1 


(4-13) 


||  6z  -  A* 


(4-14) 
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If  equation  (H-31)  is  substituted  into  equations  ( U— 23)  and  (4-2l*)s 
one  may  solve  for  the  multipliers  v  and  v.  These  are  given  by 
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V  =  -(A4*  w"1A4’)'1{2cA4>  +  A*  w-1AJ 

-  A4*  w“1A4'M“J[2cAvi  -  2CA*  w"1A4'(A*  v-1  A*)-1A*  +  M^]}  ( U-32 ) 

and 

T  T 

v  =  -M7*{2ui  -  2CA4,  w^A^A*  w'Vr1**  +  M  _}.  (U-33) 

Vw  yu 

The  matrix  M ,  .  and  the  vector  M . _  are  defined  as  follows : 

M44  =  ^  -  A*(A*  v-1A^ )_1A^  w“1]a’1'  (4-34) 

T  T  T 

M  _  -  A’4*  w_1[I  -A*(A*  w*1A4,)_1A4’  w-1]AJ.  (4-35) 

If  is  empty,  M  is  not  defined.  Similarly,  when  $  is  empty  the 

fp  T  T 

(f) A  —1  d) 

product  A9  v  A9  is  not  defined,  so  that  provisions  must  be  made  to 

properly  define  equations  (4-32)  through  (4-35)  for  all  cases.  This 

can  be  done  using  the  following  definitions.  If  $  iB  empty,  set 

T 

A<}>  =  0,  A^  e  0,  and  A4*  w  ^A4  =  I.  If  4>  is  empty,  put  £4  *  0,  A4»  “  0, 

and  =  I.  Jhese  definitions  will  reduce  equations  (4-32)  through 

(4-35)  appropriately.  The  multipliers  (4-32)  and  (4-33)  are  then 
used  to  write  the  required  change  in  the  design  variables  in  the  fol¬ 
lowing  form: 

6xS~k6X!J'  6X2’ 

where 

fix..  =  w-1[I  -  A*(A*  w^A4 )_1A4  w-1](AJ  -  A4'm7^Mit)  (4-36) 

1  1 tf/J 
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The  vectors  6x^  and  5x^  have  the  following  properties: 


fiXjWfiXg  =  0 


Av  6xr 


A4> 


Av  6x2  =  At}> 
/«*,  =  0 


Av  6Xj^  =  0 


A  6xx  >  0 


*2 


(*-37) 

(*-38) 

(*-39) 

(*-*0) 

(*-41) 

(*-*2) 

(*-*3) 


Fran  these  properties  and  from  (*-36)  it  is  seen  that  6x^  is  a 
projection  of  the  gradient  vector  on  the  planes  tangent  to  the  con¬ 
straints  and  thus  reduces  the  value  of  the  objective  function  con¬ 
sistent  with  the  constraints.  The  vector  6Xg  is  a  correction  vector 
which  directs  the  solution  towards  the  feasible  region.  Furthermore 
these  two  vectors  are  orthogonal  with  respect  to  the  weighting  matrix 
v  and  are  therefore  in  that  sense  independent. 

The  properties  of  6x1  and  6x2,  and  equation  (*-25),  can  be 
used  to  determine  the  parameter  5. 


'  ^-^jrr^T-r-T-r.-.f-i^^-.- 


E  ^ 
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6xTw6x  »  6x!fw6x..  +  $x;Fw<5x_  =  52  (4-44) 

k?  1  1  22 

Fran  this  expression  it  is  clear  that  ?  is  not  arbitrary.  In  fact,  it 
depends  on  5  which  from  (4-44)  must  satisfy 

52  >_  SXgWoXg.  (4-1*5) 

In  addition,  if  there  are  no  constraints  or  state  variables,  the 
expression  for  6x  reduces  to 


and  the  move  to  a  new  approximation  of  the  minimum  is  in  the  direction 
of  steepest  descent.  Therefore,  it  is  required  that  ;  >  0.  One  pro¬ 
cedure  that  can  be  applied  to  determine  the  size  of  step  made  in  the 
direction  of  Sx  is  to  choose  £  and  then  use  (4-44)  to  compute  c. 
Alternatively,  one  can  use  the  following  expression, 

Sx  »  -n6x1  +  5Xg 

and  choose  n  >  0  and  small.  The  problem  of  choosing  the  step  size 
will  be  considered  further  in  the  next  section  along  with  other  com¬ 
putational  considerations . 

The  constrained  steepest  descent  technique  has  been  developed 
for  problem  (4-2).  The  optimal  structural  design  problem  being  con¬ 
sidered  in  this  research  has  been  cast  in  the  form  of  problem  ( 4—1 ) 
which  is  included  in  problem  (4-2)  as  a  special  case.  All  that  is 
required  is  to  note  that,  since  the  objective  function  and  the  <j> 


constraints  are  not  functions  of  the  state  variables,  XJ  *  0  and 
^  ~  0  from  equations  (4-l6)  and  ( it— 17 ) . 

This  simplifies  the  expressions  for  AJ  and  —  equations  ( 4-26 ) 
and  (4-27).  The  remaining  expressions  are  not  changed. 

4 . 3  Computational  Considerations 

The  general  procedure  for  determining  the  optimum  of  problem 
(4-1)  or  (4-2)  is  outlined  below. 

(1)  Make  an  estimate  of  the  optimal  design  vector,  x^. 

(2)  Solve  for  the  state  variables,  z,  corresponding  to  the 
design  vector  of  the  current  iteration. 

C'  Solve  (4-16),  (4-17),  and  (4-18)  for  AJ,  1*,  and  A*. 

(4)  Determine  AJ,  A^ ,  and  A^,  from  equations  (4-26),  (4-27), 
and  (4-28). 

(?)  Choose  A$  and  A<p  and  compute  M.  .  and  M .  _  free  equations 
(4-34)  and  (4-35). 

(6)  Compute  6x1  and  6xg  from  (4-36)  and  (4-37). 

(7)  Choose  n  >  0  and  compute 

x^+1  *  xJ  -  n6x1  +  SXg. 

(8)  Check  for  convergence  and  terminate  or  go  to  (2). 

There  are  several  points  in  the  above  algorithm  that  require  further 

comment.  First  of  all,  in  the  course  of  computations  two  matrices, 

T 

.  *  -i  *. 

(A  w  A  )  and  M^,  have  to  be  inverted  and  therefore  must  be  non- 
singular.  Secondly,  the  quantities  A+  and  Atj/  must  be  chosen.  A 
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third  consideration  is  the  choice  of  the  appropriate  step  size  as 
determined  by  the  choice  of  n,  and  a  fourth  concern  is  related  to  the 
tendency  of  the  method  to  keep  a  constraint  satisfied  once  it  becomes 
identically  zero. 

The  tvo  matrices  of  concern  are  positive  seal-definite.  In  the 
computational  procedure  they  are  assumed  to  be  positive  definite. 
Experience  with  the  algorithm  has  shown  that  this  assumption  is  almost 
always  valid.  Certain  pathological  cases  can  arise,  however,  when 
this  is  not  true.  An  examination  of  equation  (4-34)  shows  that  M  ^ 
will  be  singular  when  the  columns  of  A  are  linearly  dependent.  This 
follows  from  the  fact  that  the  rank  of  a  product  of  two  matrices  can¬ 
not  exceed  the  rank  of  either  factor.  Similarly,  linear  dependence 

T 

between  columns  of  A1*'  will  make  the  product  (A^  v*A*)  singular. 

Examination  of  equations  (4-17),  (4-18),  (4-27),  and  (4-28)  shows 

(b  \i» 

that  the  matrices  A‘  and  A  are  ,:losely  related  to  the  gradients  of 
the  $  and  constraints.  For  example,  when  the  4>  constraints  are  not 
dependent  upon  the  state  variables  (as  in  the  structural  problem  con¬ 
sidered),  then  the  columns  of  the  matrix  are  just  the  gradient 
vectors  of  the  violated  or  satisfied  constraints.  If  these  gradi¬ 
ent  vectors  are  linearly  dependent  then  it  follows  di..*ectly  that 
T 

/  4  -1  4, 

(A  w  A  }  is  singular.  Furthermore,  the  assumptions  in  the  Kuhn- 
Tucker  necessary  conditions  are  not  satisfied.  The  matrix  can 

also  become  singular  if  A  is  a  square  matrix.  This  will  happen 
when  the  number  of  $  constraints  violated  is  equal  to  the  number  of 
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design  variables.  The  direct  consequence  is  that  the  product 

T  T 

tor  4  -1  4,-1  4  -1 

A  (A  w  A  )  A  w  becomes,  the  identity  matrix  making  =  0. 
Other  such  cases  could  probably  be  identified.  The  author's  experi¬ 
ence  indicates  that  these  are  the  exception,  rather  than  the  rule, 
and  resulted  from  the  introduction  of  redundant  constraints  or  a  poor 
choice  of  n  in  determining  step  size.  In  all  cases  the  cause  of  the 
singularity  was  readily  identified  and  corrected  without  altering  the 
basic  algorithm. 

Equation  (4-37),  which  is  used  to  compute  6x^,  requires  the 
quantities  A  4  and  A4.  The  vector  6Xg  is  the  correction  component  of 
the  vector  Sx  which  directs  the  search  for  the  optimum  back  towards 
the  feasible  region.  The  magnitude  of  this  correction  depends  upon 
the  values  assigned  to  the  vectors  A4  and  A$.  In  the  structural 
design  problem  the  values  A$  »  and  A$  =  were  assigned.  This 
choice  is  simple  and  effective. 

The  value  of  n  controls  both  the  direction  and  magnitude  of  the 
change  in  the  design  variable  vector  ox.  Several  methods  of  choosing 
n  were  tried.  The  most  successful  of  these  was  to  specify  a  maximum 
value  of  n  and  to  use  n  =  unless  the  constraints  were  violated 
■beyond  certain  preassigned  limits.  For  example,  for  a  maximum  allow¬ 
able  stress  of  40,000  psi,  n  *  was  used  unless  the  stresses 
obtained  using  this  value  exceeded  44,000  psi.  In  other  words  a  ten 
per  cent  violation  of  the  constraints  was  allowed.  This  sped  up  con¬ 
vergence  by  allowing  bigger  steps  to  be  taken  without  excessively 
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violating  the  constraints.  If  an  unacceptable  violation  of  the  con¬ 
straints  occurred,  the  current  vaj.ue  f  n  was  reduced  by  ten  per  cent 
and  the  constraints  were  checked  again.  This  process  was  continued 
until  an  acceptable  value  of  n  was  obtained  or  until  n  became  zero. 

If  this  occurred,  the  value  of  n  was  set  no  10”^  and  computations  were 
allowed  to  proceed.  A  zero  value  of  n  is  not  allowed,  since  at  oti*er 
points  in  the  algorithm  division  by  n  is  required. 

When  an  estimate  of  the  solution  is  such  that  seme  of  the  con¬ 
straints  are  exactly  satisfied,  the  method  demands  that  they  continue 
to  be  satisfied.  It  may  be  advantageous  if  the  search  for  the  mini¬ 
mum  is  allowed  to  leave  such  a  constraint.  Theorem  4.1  Btates  that 
the  multipliers,  u  and  v,  which  exist  at  a  solution,  must  be  positive 
or  zero  tut  that  they  cannot  be  negative.  Therefore,  at  each  suc¬ 
cessive  iteration  of  the  optimisation  process  the  multiplier  vectors 
y  and  v  are  computed.  If  either  of  these  has  components  which  are 
negative,  the  constraints  corresponding  to  these  components  are 
removed  from  the  appropriate  constraint,  set  £  or  This  procedure 
proved  effective  when  the  method  was  applied  to  the  optimal  structural 
design  problem. 


4,4  Convergence  Cr iter- ion 

A  convergence  criterion  must  be  established  for  Step  (8)  in  the 
algorithm  of  Section  4.3.  It  will  be  shown  that,  if  the  sequence  of 
solutions  of  che  linearized  problem  converges  to  a  ’ical  solution  of 
the  nonlinear  problem,  then  <Sx^  *  0.  The  result  is  obtained  by  the 
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direct  application  of  the  Kuhn -Tucker  necessary  conditions  to  both  the 
nonlinear  and  the  linearized  problems  and  comparison  of  the  resulting 
equations. 

Theorem  4.1  states  that  at  the  solution  of  the  nonlinear  problem 
(4-2)  there  exist  multipliers 


^  Q,  X  *  •  •  •  g  £* 

and 

i  **  Xf  •••§  it 

such  that  for 

H  =  f(z,x)  +  XTh(z,x)  +  yT^(z,x)  +  vT\j»(z) 

then 

3H  3f  .T  3h  T  3$  _ 

T**  *  r~  +  X  t—  +  u  r**  **  0 

3x  3x  3x  3x 

(4-46) 

3H  3f  T  3h  ,  13*  T3t  _ 

3z  3z  3z  3z  3z 

(4-47) 

Mi  =0  if  ^(z.x)  <  0 

and 

=  0  if  li..  (*)  <  0. 

Define  *  and  $ 

such  that 

?  ~  Uji  ^U,*)  >  0} 

(4-48) 

$  a  {*  •  V^(z)  5^  0} 

(4-49) 

and  y  and  v  which  contain  only  the  components  of  y  and  v  respectively, 
corresponding  to  the  constraints  in  $  and  Using  this  information, 
equations  (4-46)  and  (i»-47)  can  be  written  in  the  following  form. 
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g.X*£.?g.. 

£*X*£.  ?&*?£... 


(4-50) 

(4-51) 


Using  equations  (4-l6)  through  (4-13) ,  equation  (4-51)  may  he 
written 

fc*1  a  +  $*/  a  .  (X-52) 


3h 


The  matrix  Is  nonsingular,  therefore 

AJ  +  X  ♦  1*5  ♦  X^v  -  0.  (4-53) 

Equation  (4-53)  may  be  used  to  eliminate  X  from  equation  (4-5?.)  to 


obtain 

g  -  <ajT  ♦  ;/  *  ;VT)  g  ♦  51  f§  .  o,  (W5M 

Collecting  coefficients  on  y  and  v,  and  using  equations  (4-2 6)  through 
(4-28)  yields, 

A*y  »  -(AJ  «■  A*u),  (4-55) 

*T  -1 

If  equation  (4-55)  i»  preavltiplied  by  Af  v  ,  the  following  expression 
for  y  is  obtained 


V  .  -(a*  v^A+rV  v"1(AJ  +  A*C)  (4-56) 

-l  a 

It  should  be  noted  that  it  has  been  assumed  that  (A*  v  A*)  is  non- 
singular.  If  equation  (4-56)  is  substituted  into  (4-55), 

[i  -  A*(A*  v”1A^)“'1’A*  w"^3[aJ  +  A*y]  ■  0  (4-57) 

T 

A*  J 

is  obtained.  Premult iplying  by  Av  v  and  using  equations  (4-34)  and 
(4-35),  equation  (4-57)  becomes 


» 


Therefore  by  equation  (U-36),  at  the  solution  of  problem  (4-2),  it  is 
necessary  that  ta^  *  0.  Furthermore  in  the  limit  +  ■  0  and  *  ■  0 
so  that  A#  «  0  and  hi  *  0  and  by  equation  (fc-37)»  tag  *  0.  Since 

ta  «  -nta^  +  tag,  (H*-60) 

ta  must  also  approach  zero  if  the  procedure  converges.  Satis  fiction 
of  this  criterion  is  used  in  the  algorithm  given  in  thi  next  section. 


51 


4.5  CSDS  Algorithm 

The  following  algorithm  is  an  expansion  of  the  procedure  defined 
in  Section  4.3  and  incorporates  the  additional  features  discussed  in 
Sections  4.3  and  4.4. 

(1)  Estimate  the  optimal  design  vector,  rS^ . 

(2)  Determine  the  state  variables  z  corresponding  to  the 
design  variables  of  the  current  iteration,  x^. 

(3)  Solve  (4-l6),  (4-17),  and  (4-18)  for  AJ,  X*,  and  A*. 

(4)  Determine  $  and  iji. 

(5)  Determine  A1^,  A^,  and  A^  from  equations  (4-26),  (4-27), 


and  (4-28). 

v.6)  Choose  Aip  =  and 

(7)  Compute  and  from  equations  (4-34)  and  (4-35). 

(8)  Determine  the  multipliers  y  and  v  from  (4-32)  and  (4-33). 

(9)  If  any  of  the  components  of  y  and  v  are  negative,  remove 

«*  ** 

the  corresponding  constraints  from  the  sets  $  and  ♦  and 
return  to  Step  (5).  If  all  components  of  y  and  v  are  non¬ 
negative,  continue. 

(10)  Compute  6x^  and  &x2  from  (4-36)  and  (4-37). 

(11)  Choose  n  =  nmPT. 

(12)  If  n  ^  0,  set  n  =  10~^  and  go  to  Step  (15)' 

(13)  Compute 

~  x^  -  +  dx2. 
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(lU)  Compute  the  values  of  the  constraints  for  x^+\  If  the 
constraints  have  been  violated  excessively,  reduce  n  and 
return  to  Step  (ll) ,  otherwise  continue. 

(15)  If  <5x^  <ex»  i»l,  ...»  n  terminate.  Otherwise  return 
to  Step  (2). 

The  application  of  this  algorithm  to  the  optimal  structural 
design  problem,  and  the  ensuing  results,  sore  discussed  in  Chapter  V. 


CHAPTER  V 


APPLICATIONS 


5.1  General  Discussion 

The  optimization  methods  presented  in  Chapter  III  and  Chapter  IV 
wer^  applied  to  the  design  of  several  two  and  three  member  frames  sub¬ 
jected  to. a.  variety  of  loads.  In  order  to  aid  discussion  the  frame 
members  and  loads  are  numbered  as  shown  in  Figure  5.1  and  Figure  5.2. 
All  dimensions  and  limits  on  dimensions  are  given  in  inches,  the  loads 
are  in  pounds  and  the  stresses  are  in  pounds  per  square  inch.  In 
Section  5.2  and  5.3  a  general  descr  ption  of  the  frames  to  be  opti¬ 
mized  is  given.  Sections  5.^  and  5.5  present  a  brief  discussion  of 
the  programs  used  and  some  aspects  of  their  application.  The  result¬ 
ing  designs  are  presented  and  discussed  in  Sections  5.6  and  5.7. 

5.2  The  Two  Member  Frame 

The  two  member  frames  vhich  were  optimized  are  shown  _n  Tables 
5.1  through  ''.IG.  In  each  case  there  were  sis  design  'mr^abl^s  to 
be  determined,  ten  state  variables,  and  seventeen  constraints  to  be 
satisfied.  The  design  variables  sure  the  vail  thickness,  t,  width,  b, 
sind  height,  h,  (see  Figure  2.2)  of  each  member.  These  variables  must 
satisfy  the  following  inequalities  for  the  i*1*1  member: 


Member  1 


—  Member  2 


FIGURE  5*1  Numbering  the  T to  Member  Frame 


FIGURE  5.2  Numbering  the  Three  Member  Frame 


0.1  <_  ti  <_  1.0 

2.5  1^1  10*° 
2.5  <  h.  <  10.0. 

mm  1  «M» 


(5-1) 


The  state  variables  are  the  three  displacement  components  of  the 
joint,  the  torsional  stress  in  each  member,  and  the  bending  stresses 
at  the  five  critical  points  of  the  structure.  The  critical  points 
are  the  ends  of  each  member  and  under  the  applied  load.  The  stresses 

JL-L. 

must  satisfy  the  following  inequality  at  the  i  critical  point: 


"  *0,000  1  0 


(5-2) 


where  is  the  normal  stress  due  to  bending  and  is  the  shear 
stress  due  to  torsion  at  the  point. 

The  sequential  unconstrained  minimization  technique  and  the 
method  of  constrained  3teepest  descent  with  state  equations  were  both 
applied  to  the  two  member  frames.  The  comparative  effectiveness  of 
these  two  approaches  iB  discussed  in  Section  5.6» 


5.3  The  Three  Member  Frame 

The  three  member  frames  optimized  are  shown  in  Tables  5.11 
through  5«1T.  For  each  frame  nine  design  variables  were  chosen  sub¬ 
ject  to  eighteen  design  variable  constraints.  There  are  eighteen 
state  variables  and  twelve  state  variable  constraints.  The  design 
variables  are  the  same  as  those  for  the  two  member  frame.  The  state 
variables  consist  of  six  joint  displacements  and  twelve  stresses. 

The  stresses  are  th*  torsional  stress  in  each  member  plus  three 
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critical  bending  stresses  in  each  member  occurring  at  the  ends  or 

under  the  applied  load.  For  the  frames  in  Tables  5.11  through  5*15* 

till 

the  stresses  at  the  i6  critical  point  were  required  to  satisfy  the 
condition  (5-2),  and  the  design  variables  were  required  to  satisfy 
(5-1). 

The  dimensions  and  loads  shown  in  Tables  5.16  and  5.17  were 
chosen  to  approximate  a  frame  that  the  author  encountered  in  an 
industrial  application.  In  both  cases  the  stresses  at  a  critical 
point  are  required  to  satisfy  the  following  condition: 

^oj  +  3t|  -  20,000  <.  0.  (5-3) 

The  design  variables  for  the  frame  in  Table  5.15  satisfy 
.109  <  t±  <  1.0 

2.0<bi<  iO.O  (5-h) 

6.0  -  hi  1  10*°  for  i  =  1»  3, 

and 

.25  <.t2  <1.0 

3.°  1  -  10,0  (5-5) 

3.0  <  h„  <  10.0. 

The  design  variable  constraints  which  the  frame  in  Table  5.16  satis¬ 
fied  are, 

.109  <,  <1.0 

2.0  <  b±  ^  5.0 
6.0  <h±<  12.0 

for  members  1  end  3,  and  (5-5)  for  member  2, 


(5-6) 
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Because  of  its  superior  perforrance  on  the  two  member  frames, 
omy  the  method  of  constrained  steepest  descent  with  state  equations 
was  used  in  the  optimum  design  of  the  three  member  frames. 

5 .4  The  SUMT  Program 

The  sequential  unconstrained  minimization  technique  was  pro¬ 
grammed  using  the  algorithms  of  Sections  3.1  and  3.2.  Acceleration 
by  extrapolation  was  included  and  proved  to  be  effective  both  in 
decreasing  the  number  of  function  evaluations  and  in  giving  better 
convergence  to  the  optimum,  but  only  when  the  calculations  were  done 
in  double  precision.  Fiacco  and  McCormick  elude  to  this  when  they 
discuss  the  significant  effect  of  round-off  error  in  the  extrapola¬ 
tion  procedures. 

Certain  safeguards  must  he  introduced  into  the  program  to  insure 
that  the  search  for  a  minimum  always  takes  place  within  the  feasible 
region  as  required  by  the  development  in  Chapter  III.  In  the 
Fibonacci  search  increasingly  larger  steps  are  taken  until  the  mini¬ 
mum  is  bracketed.  When  the  minimum  is  near  the  boundary  of  the  fea¬ 
sible  region,  it  is  possible  for  a  step  to  violate  a  nearby  con¬ 
straint.  Therefore  any  time  a  new  step  is  taken  in  the  optimization 
process  a  check  must  be  made  to  see  if  a  constraint  has  been  violated. 
If  violation  occurs,  then  appropriate  measures  must  be  taken  to  return 
the  search  to  the  feasible  region.  These  procedures  axe  ad  hoc  and 
significantly  reduce  the  efficiency  of  the  unconstrained  minimization. 


i 


Three  convergence  criteria,  c  ft,  and  e..  *  must  be  chosen  for 

sx  p  j» 

SIMT.  Of  these  should  have  the  ranaileat  value,  since  the  greatest 

accuracy  is  required  in  the  minimization  of  a  function  along  a  line. 

e  cannot  be  less  than  the  values  of  either  e  or  eT .  The  choice  of 
sx  p  L 

these  parameters  significantly  affects  the  convergence  of  the  method 
and  in  almost  all  cases  experimentation  was  required  before  a  satis¬ 
factory  set  was  obtained. 


5.5  CSDS  Program 

The  constrained  steepest  descent  program  was  written  using  the 
algorithm  of  Section  U.5.  All  computations  are  done  in  double  pre¬ 
cision  and  the  program  is  written  so  that  the  required  derivatives 
can  be  computed  exactly  or  from  finite  difference  approximations. 
Derivatives  of  the  state  equations,  constraint  functionsB  and  objec¬ 
tive  functions  are  required  by  the  algorithm.  The  derivatives  of  the 
state  equations  with  respect  to  the  design  variables  are  the  only 
derivatives  approximated  since  the  remaining  functions  are  simple  and 
the  derivatives  can  be  computed  exactly  without  undue  effort.  Each 
frame  shown  in  the  tables  wan  optimized  using  both  exact  and  approxi¬ 
mate  derivatives,  A  simple  forward  difference  scheme  was  used  ini¬ 
tially  to  obtain  the  approximate  derivatives.  Some  difficulty  was 
experienced  but  results  were  improved  when  the  interval  of  the  approx¬ 
imation  was  reduced.  Still,  for  seme  three  n  amber  frames,  the  method 
did  not  converge  properly  when  the  derivatives  were  approximated. 
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If 


This  problem  was  eliminated  at  the  expense  of  more  function  evalu¬ 
ations  by  using  a  central  difference  approximation. 

Two  parameters,  and  ,  must  be  chosen  when  the  program  is 
used,  n  is  the  ma-^-tranm  step  size  allowed  and  e  is  the  convergence 

mav  r  X 

Ji 

criterion.  A  value  of  e  =  .1x10  was  chosen.  This  is  somewhat 

x 

more  stringent  than  the  value  of  used  is  SUMT;  but  if  the  method 

converges  at  all,  it  has  no  difficulty  satisfying  this  criterion. 

Some  minimal  experimenting  may  he  necessary  to  choose  n  .  For  large 

values  of  n  ,  fewer  iterations  are  required  if  the  method  converges, 
max* 

hut  sometimes  situations  may  be  created  which  the  algorithm  cannot 
correct.  These  sire  immediately  apparent  and  easily  eliminated  by 
reducing  the  value  of 

5.6  Application  to  Two  Member  Frames 

Three  different  design  problems  are  shown  in  Tables  5.1  througn 
5.10.  These  designs  were  optimized  using  both  the  sequential  uncon¬ 
strained  minimization  technique  and,  the  method  of  constrained  steepest 
descent  with  state  equations.  The  frame  shown  in  Tables  5.3  through 
5.4  has  members  of  different  lengths  with  the  loud  applied  at  the 
joint.  The  computation  time  for  the  solution  using  SUMT  {Sable  5.1) 
is  considerably  greater  than  the  time  required  using  CSDS  (Tables  5*2 
through  5.3),  even  though  the  convergence  criterion  e  for  SUMT  is 
considerably  less  restrictive  than  the  criterion  used  in  CSDS.  In 
SUMT  the  state  of  the  structure  is  computed  for  each  function  evalu¬ 
ation.  In  CSDS  the  state  Is  computed  about  ten  times  per  iteration. 


:!! 
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Using  this  value  the  number  of  functional  evaluations  required  for 
the  design  of  Table  5*3  is  about  140  compared  to  the  8,633  ftvalu- 
aticns  required  for  SUMT.  The  results  of  Tables  5.1  through  5.3 
chcv  thet  different  starting  values  converge  to  the  same  minimum. 
However,  when  the  starting  values  of  Table  5.1  were  used  in  the  CSDS 
method,  another  local  minimum  was  obtained.  This  is  not  surprising 
in  view  of  the  nonconvex  nature  of  the  programming  problem. 

Another  frame  was  optimized  using  SUMT  (Tables  5.5  through  5.6) 
and  CSDS  (Table  5.7).  The  large  difference  in  computing  time  for  the 
two  methods  is  again  apparent.  Starting  from  the  same  values,  the 
two  methods  converged  to  two  different  local  minima.  Several 
attempts  were  made  to  start  in  the  neighborhood  of  the  m-jn-tmnm  deter¬ 
mined  by  SUMT  using  CSDS,  but  the  latter  technique  always  converged  to 
the  symmetric  results  of  Table  5.7.  An  additional  trial  (Table  5.6) 
was  made  using  SUMT  with  a  starting  value  in  the  neighborhood  of  the 
symmetric  solution.  The  results  are  close  to  those  of  Table  5.7. 

The  volume  is  smallest  for  the  symmetrical  design. 

The  third  two  member  frame  that  was  optimized  is  shown  in  Tables 
5.8  through  5.10.  Two  trials  using  SUMT  and  one  using  CSDS  were  made. 
The  resulting  volumes  do  not  differ  greatly,  but  the  design  variables 
show  significant  differences.  About  840  function  evaluations 
(assuming  ten  per  iteration)  are  required  when  CSDS  is  used  compared 
to  11,610  required  by  SUMT  (Table  5.9) c  The  amount  of  computing  time 
required  shows  a  similar  advantage  for  CSDS. 


The  results  show  that  the  method  of  constrained  steepest  descent 
with  state  equations  has  significant  advantages  over  the  sequential 
unconstrained  minimisation  technique.  Both  the  computation  time  and 
the  number  of  times  that  the  state  of  the  structure  must  be  evaluated 
are  considerably  less  than  those  required  for  SUMT.  Furthermore,  the 
same  results  are  obtained  with  or  without  the  exact  calculation  of 
derivatives. 


P  »  10,000  lb 


of  =  39,945  psi.  at  1 


Initial 

Optimu 

Values*  in* 

Values*  ^n* 

csx  ■  ■ 1x10-2 

*1 

0.9 

0.1010 

ex  »  .  Ixl0_1+ 

bl 

9.0 

a.  5070 

e1  x  .lxlO-5 

hl 

9.0 

2.5070 

*2 

0.9 

0.1012 

b2 

9.0 

9.1*234 

h2 

9.0 

9.9827 

Volume,  in? 

4,374.00 

291.59 

Number  of  iterations 

6 

Computing  time,  sec. 

370 

Number  of  function  evaluations 

8,633 

TABLE  5*1  Optimum  Design  of  a  Two  Member  Frame  Using  SIWT 


p  «  10,000  lb 


100"  - 


40,000  psi.  at  1 


n  a  .1 
max  . 

cx  *  .1x10 


Initial 

Values, 

in. 


0.10 

2.50 

2.50 

0.11 

10.00 

10.00 


Volume,  in? 


Humber  of  iterations 


Computing  tine,  sec. 


TABLE  5.2  Optimum  Design  of  a  Two  Member  Frame  Using  CSDS 


Optimum  Values ,  in. 

MU  i 1  M  l  mu  H  1  —  1 

Exact 

Derivatives 

0.10000 

0.10000 

2.50000 

2.50000 

2.50000 

2.50000 

0.10000 

0.10000 

9.52614 

0.52614 

10.00000 

10.00000 

289.26 

289.26 

9 

9 

20 

20 

4 


TVax  * 

-4 

e  »  .1x10 

X 

*1 

\ 

hl 

*8 

*2 

h2 

Initial 

Values, 

in. 

Optimum  Vi 
Obtained! 

ilue?,  in. 
Jsins: 

Sxact 

Derivatives 

Approximate 

Derivatives 

O.UTfcl 

2.5246 

2.5217 

0.22061 

3.2585 

3.2544 

0.1C000 

2.50000 

50000 

0.10000 

9.52614 

10.00000 

0.10000 

2.50000 

2.50000 

0.10000 

9.52614 

10.00000 

Volume,  in? 

246.93 

289.26 

289,26 

Humber  of  iterations 

14 

14 

Computing  time,  sec. 

14 

14 

TABLE  5-3  Optimum  Design  of  a  Two  Member  Frame  Using  CSDS 


ico\ 

40,000  psi,  at  1,2 


*  .1x10 


P  ■  10,000  lb. 

I  . 


Initial 

ValueB 


Optimum  ValueB,  in. 
Obtained  Using 


rivatives 


0.11373 

10.00000 

10.00000 

0.23599 

10.00000 

2.5000C 


0.11373 

10.00000 

10.00000 

0.23598 

10.00000 

2.50000 


Volume,  in? 

i  4,374.00 

733.59 

Humber  of  iterations 

19 

Computing  time,  sec. 

18 

TABLE  5.4  Optimum  Design  of  a  Two  Member  Frame  Using  CSDS 


P  c  10*000  lb 


•  Initial 

Optimum  Values,  in. 
Obtained  Using 

Values  * 
in. 

Exact 

Derivatives 

Approximate 

Derivatives 

^max  = 

tl 

0.9 

0.10000 

0.10000 

=  .1x10“** 

bl 

9.0 

9.5H738 

9.5^738 

hl 

9.0 

10.00000 

10.00000 

*2 

0.9 

0.10000 

0.10C00 

*2 

9.0 

9.5UT38 

9.51*738 

• 

h2 

9.0 

10.00000 

10.00000 

Volume,  la? 

5,932.00 

773.89 

773.89 

Number  of  iterations 

15 

15 

Computing  time,  sec. 

22 

22 

TABLE  5.7  Optimum  Design  of  a  Two  Member  Fl*«ar.e  Using  CSDS 
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Initial 

Optimum 

Values,  in 

Values,  in. 

e  =  .  lxlQ**1 
sx 

tl 

0.2 

0.100 

Ex  »  . Ixl0~3 

bl 

5-0 

U.319 

=  .lxlO-1* 

hl 

9.0 

9.253 

CM 

-P 

0.2 

0.100 

b2 

3.0 

2.502 

h2 

9.0 

9.075 

Volume,  in? 

776.00 

362.07 

Number  of  iterations 

9 

Computing  time,  sec. 

310 

Number  of  function  evaluatio.  s 

11,251 

TABLE  5*8  Optinnsn  Design  of  a  Two  M saber  Frame  Using  SUMT 


f 


Initial 


Optimum  Values,  in 
Obtained  Using 


TV  j  VVW  V*  V 


! 

Values, 

in. 

1 

333531 

n  -  .1 

max 

*1 

0.9 

0.10000 

0.10000 

e  =  . lxlO-^ 

X 

bl 

9.0 

h.  1*0827 

fc.b0827 

hl 

9.0 

8.58807 

8.58807 

*2 

0.9 

0.10000 

0.10000 

*2 

9.0 

2.50000 

2.50C00 

1  _ 

h2 

9-0 

8.07038 

8.07038 

Volume,  in? 

If,  37^  .00 

359.63 

359.63 

Humber  of  iterations 

8U 

8k 

Computing  time,  oec. 

72 

_IL_1 

TABLE  5.10  Optimrai  Design  of  a  Two  Member  Frame  Using  CSDS 


The  results  of  Section  5.6  show  the  advantages  of  the  method  of 
constrained  steepest  descent  with  state  equations  over  the  sequential 
unconstrained  minimization  technique  when  optimization  without  com¬ 
puting  exact  derivatives  is  required.  Furthermore,  as  the  number  of 
design  variables  is  increased,  SUMT  will  become  increasingly  more 
difficult  to  apply  because  of  the  nature  of  the  algorithms  used  to 
determine  the  unconstrained  minima.  Therefore,  the  optimization  of 
the  three  member  frames  was  carried  out  by  using  only  CSDS, 

A  variety  of  symmetrical  and  onsymetrical  load  sets  were 
applied  to  several  different  three  member  frames.  The  minimum  weight 
designs  obtained  are  shown  in  Tables  5*H  through  5.17  *  The  agree¬ 
ment  between  designs  obtained  using  exact  and  approximate  derivatives 
is  particularly  noteworthy.  This  was  improved  during  the  course  of 
the  applications  by  introducing  a  central  difference  formula  to 
approximate  the  derivatives  of  the  3tate  equations.  The  improvement 
that  resulted  may  be  seen  by  comparing  the  results  of  Tables  5.13 
and  5.14.  The  designs  of  Tebles  5.15  through  5.17  also  show  the 
value  of  this  modification,  since  the  results  for  these  frames  using 
exact  and  approximate  derivatives  agree  exactly. 

The  application  in  Table  5.12  illustrates  another  very  useful 
feature  of  CSDS.  In  this  case  the  method  did  not  converge  within  the 

maximum  allowable  iterations,  but  the  value  of  ox  had  decreased  and 

-2 

was  of  the  order  10  .  Thus,  the  necessary  condition  t)At  ox  =  0 
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(Section  4.U)  is  approximately  satisfied,  indicating  that  the  results 
are  probably  in  the  neighborhood  of  a  minimum  weight  design.  Conse¬ 
quently,  further  iterations  are  not  absolutely  necessary.  Another 
feature  of  CSDS  shown  by  the  results  is  that  constraints  may  be  satis¬ 
fied  exactly.  In  SUMT  (see  Section  5.6)  constraints  are  only  approx¬ 
imately  satisfied. 

An  unsymmetrical  frame  (Table  5.15)  required  several  attempts  to 
obtain  the  design  shown.  This  is  not  uncommon  in  that  'jorae  experi¬ 


mentation  is  almost  always  required  to  choose  ex  and  nrnflT  properly. 
Furthermore  an  unfavorable  starting  point  may  be  inadvertently  chosen. 

The  designs  presented  in  Tables  5<l6  and  5.17  illustrate  the 
dependency  of  the  results  both  on  the  configuration  of  the  basic 
frame  and  the  assignment  of  the  design  variable  constraints  (see 


Section  5.3). 


TABU5  5*11  Optimum  Design  of  a  Three  Member  Frame  Using  CSDS 


% 


J 

,3* 

#  $ 


.  I  * 

\ 

k 


IfC 


100" 

Oj.  *  Uo,000  psi.  at  1,2 


e  =  .1x10 
x 


P  «  10,000  lb. 


Initial 

Values, 

in. 


Volume,  in?  8,7U8.Q0  1,656.91  1,656, 

Number  of  iterations  lU  lU 

Computing  time,  sec.  US  U8 

TABLE  5.lU  Optimuu  Design  of  a  Three  Member  Frame  CSDS 


Optimum  Values,  in. 
Obtained  Using 

Exact 

Derivatives 

Approximate 

Derivatives 

0.19908 

0.19908 

10.00000 

10.00000 

10.00000 

10.00000 

0.10000 

0.10000 

2.50000 

2.50000 

2.50000 

2.50000 

0.19908 

0.19908 

10.00000 

10.00000 

10.00000 

10.00000 

1,656.91 

1,656.91 

lU 

lU 

US  | 

U8 

Initial 

Optimum  Values,  in. 
Obtained  Using 

Values t 
in. 

Exact 

Derivatives 

(Vpproximate 

Derivatives 

n  *  . lxlO-1 

max 

’‘1 

0.9 

0.10000 

0.10000 

e  n  .lxlO-i+ 

X 

bl 

9.0 

9.53870 

9.53870 

hl 

9.0 

10.00000 

10.00000 

*2 

0.1 

0.10000 

0.10000 

b2 

2.5 

2.50000 

2.50000 

h2 

2.5 

2.50000 

2.50000 

*3 

0.1 

o.xoooo 

0.10000 

b3 

3.0 

2.50000 

2.50000 

h3 

3.0 

2.50000 

2.50000 

Vplume,  in? 

3,051.00 

506.77 

506.77 

Number  of  iterations 

32 

32 

Computing  time,  sec. 

80 

80 

TABLE  5,15  Optimum  Design  of  a  Three  Member  Frame  Using  CSDS 


P  «  564  lb 


Initial 

Optimum  Values,  in. 
Obtained  Using 

Values . 
in. 

Exact 

derivatives 

Approximate 

Derivatives 

n  =  ,2xl0-1 

max 

t1 

0.9 

0 . 1339-s 

0.13394 

•  m-1* 

e  =  •  .LxlO 

A. 

bl 

9.0 

10.00000 

10.00000 

hl 

9.0 

10.00000 

10.00000 

%2 

0.9 

1.00000 

1.00000 

b2 

9.0 

3.00000 

3.00000 

-12 

9.0 

3.00000 

3.00000 

*3 

0.9 

0,13394 

0.13394 

b3 

9.0 

10.00000 

10.00000 

• 

h3 

9.0 

10.00000 

10.00000 

Volume ,  in? 

5,471.00 

1,088.05 

1,088.05 

Number  of  iterations 

35 

35 

Computing  tine,  sec. 

77 

77 

TABLE  5.16  Optimum  Design  of  a  Three  Member  Frame  Using  CSDS 


f 


Initial 

Optimum  Values,  JJCw 
Obtained  Using 

■ 

Valves t 
in.  . 

Exact 

Derivatives 

Approximate 

Derivatives 

n  =  .2xl0_1 

max 

H 

0.9 

0.20687 

0.20687 

ex  **  ,1x10 

Q 

4,9 

5.00000 

5.00000  ' 

Wm 

11.0 

12.00000 

12.00000 

*2 

1.0 

0.25000 

0.25000 

*2 

10.0 

3.00000 

3.00000 

h2 

9.0 

3.00000 

3.00000 

t3 

0.9 

0.20687 

0.20687 

b3 

4.9 

5.00000 

5.00000 

h3 

11.0 

12.00000 

12.00000 

Volume,  in? 

5,067.81 

l.lUl.52 

1,141.52 

Humber  of  iterations 

28 

28 

Computing  time,  sec. 

70 

70 

TABLE  5.17  Optimum  Design  of  a  Three  Member  Fra ae  Using  CSDS 
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CHATTER  VI 
CONCLUSIONS 

The  results  of  Chapter  V  shov  that  a  class  of  spatial  structures 
can  be  optimally  designed  by  the  method  of  constrained  steepest 
descent  with  state  equations,.  In  particular,  this  method  appears  to 
be  superior  to  the  sequential  unconstrained  minimization  technique 
when  necessary  derivatives  cannot  be  computed  exactly.  The  class  of 
problems  which  can  be  solved  by  the  method  is  broad  and  includes 
those  struct  lores  that  can  he  analyzed  by  matrix  displacement  and 
finite  element  methods. 

In  large  measure,  the  effectivenecs  of  C3DS  results  from  the 
fact  that  its  formulation  so  closely  matches  the  essential  features 
of  the  original  design  problem.  This  is  a  more  natural  approach  to 
the  problem  and  leaves  the  associated  functions  in  their  most  simple  . 

j 

form.  Other  programming  methods  tend  to  complicate  the  x'ormulation.  ] 

In  SUMT  the  surface  created  by  adding  a  penalty  term  to  the  objective  i 

function  exhibits  characteristics  that  even  the  most  robust  uncon¬ 
strained  minimisation  technique,;  have  difficulty  coping  with.  In  ; 

other  approaches  >.  state  equations  are  used  to  write  all  of  the  con¬ 
straints  as  functions  of  the  design  variables;  aa  a  consequence,  com¬ 
plicated  nonlinear  function1;  are  obtained.  Another  method  defines 
the  set  of  optimization  variables  to  consist  of  the  design  variables 


t 
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plus  the  state  variables «  For  the  structures  considered  in  this 
research,  minimization  problems  of  dimension  sixteen  and  twenty- 
seven  would  result.  Such  problems  are  formidable  particularly  when 
the  minimization  must  be  accomplished  without  calculating  derivatives. 

The  method  of  constrained  steepest  dejcent  with  state  equations 
avoids  these  pitfalls  because  the  state  equations  and  constraints 
involving  the  state  variables  are  included  directly  in  the  formula¬ 
tion.  Consequently,  a  great  deal  of  flexibility  is  introduced  into 
the  solution  since  the  functions  are  retained  in  their  simplest  form. 
Therefore,  derivatives  are  approximated  only  when  it  is  absolutely 
necessary  since,  in  the  case  of  simple  functions,  derivatives  nay  be 
computed  exactly  without  .indue  effort. 

CSDS  is  an  effective  solution  to  the  optima*  structural  design 
problem  for  two  additional  reasons.  First  of  all,  much  o’  the  "art" 
is  removed ,  since  only  two  parameters  must  be  chosen  —  and 
In  SUMT  there  are  five  — e  ,  c  ,  e. ,  r1s  and  C.  Therefore,  the 

EX’'  p*  ij*  l3 

amount  of  numerical  expel’ iaentat ion  that  must  be  done  for  a  given 
problem  is  greatly  reduced.  Secondly,  the  procedure  has  been  written 
so  that  existing  structural  analysis  techniques  can  be  used  to  full 
advantage.  Consequently,  currently  available  algorithms  based  on 
matrix  displacement  and  finite  element  techniques  may  be  used  in  the 
optimization  procedure  with  minimal  effort  for  adaptation. 

In  this  research,  a  basic  configuration  of  the  structure  was 
first  chosen  and  then  the  dimensions  of  individual  members  were 


determined  using  methods  of  optimal  design.  There  are  sane  inter¬ 
esting  variations  of  this  problem  which  may  stimulate  further  inves¬ 
tigation.  These  variations  involve  both  the  structure  and  the  loading 

Structures  of  greater  complexity  can  be  generated  by  increasing 
the  lumber  of  the  members  or  increasing  the  nunber  of  design  vari¬ 
ables  per  member.  In  the  first  cast*  the  methods  developed  in  this 
research  apply  directly,  since  the  basic  form  of  the  functions 
involved  is  not  changed.  The  incorporation  of  existing  structural 
analysis  algorithms  would  be  a  natural  way  of  coping  with  the 
increased  size  of  the  problem .  An  example  of  the  second  case  is  the 
introduction  of  linear  variations  of  the  height  and  width  of  the  mem¬ 
bers.  A  problem  of  this  type  will  require  an  effective  moans  of 
determining  the  points  of  maximum  stress  in  the  structure  and  more 
extensive  numerical  work  in  the  structural  analysis. 

The  optimum  design  problem  can  be  formulated  to  include  the 
basic  configuration  of  the  structure  by  adding  the  lengths  of  the 
members  as  design  variables t  This  will  increase  the  complaxity  of 
the  functions  involved.  For  example,  the  objective  function  will  be 
a  cubic  in  the  design  variables. 

In  practice,  spatial,  structures  may  be  subject  to  a  variety  of 
randomly  applied  loads.  In  addition,  the  state  variable  constraints 
may  change  for  different  loadings.  An  optimal  design  procedure  for 
spatial  structures  with  multiple  load  sets  and  multiple  constraint 
sets  should  be  developed  for  this  type  of  problem. 
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The  method  of  constrained  steepest  descent  with  state  equations 
may  have  to  be  extended  or  modified  to  effectively  solve  seme  of 
these  problems-  For  example,  steepest  descent  methods  converge 
slowly  when  the  objective  function  has  a  long  valley.  Several 
methods  have  been  develcped  to  deal  with  problems  of  this  type  in 
unconstrained  minimization  and  perhaps  modifications  similar  to  these 
can  be  introduced  into  CSDS.  The  method  is  certainly  not  limited  to 
structural  design,  but  may  be  applied  to  a  wide  variety  of  finite 
dimensional  design  problems.  The  magnitude  and  complexity  of  the 
problems  that  this  method  can  successfully  solve  is  not  yet  known 
and  should  be  investigated  further. 
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AFPENDIX 

STRUCTURAL  AHALYSIS 


A.l  The  Basic  Equations 

The  structural  analysis  used  in  this  research  foil ova  a  method 
discussed  by  Langhaar  [11],  T’-e  analysis  of  a  two  member  frame  will 
he  done  in  detail  to  illustrate  the  application  of  the  method.  The 
equations  and  associated  matrices  required  for  the  analysis  of  the 
three  member  frame  will  then  be  shown. 

Consider  the  two  member  frame  in  Figure  A.l.  The  horizontal 
frame  is  composed  of  hollow  rectangular  members  which  are  Joined 
perpendicular  to  each  other.  The  Joint  is  rigid  and  transmits  shear 
force  as  well  as  bending  and  twisting  moments#  The  load,  P,  is 
applied  normal  to  the  frame  at  an  arbitrary  point  along  member  (1,2). 
This  point  is  denoted  in  Figure  A.l  as  point  U.  The  members  are 
clamped  at  points  1  and  ?.  The  reactions  at  these  points  are  showa 
on  the  free  body  diagram  in  Figure  A, 2.  The  deformation  of  the  frame 
is  completely  described  by  the  displacements  at  joint  2.  These  dis¬ 
placements  are  the  vertical  deflection  of  point  2,  and  the  rota¬ 
tions  of  the  joint  about  the  axes  of  the  two  members. 

The  behavior  or  state  of  the  frame  is  specified  by  the  value 
of  the  maximum  failure  stress  occurring  in  the  frame.  This  failure 
stress  is  calculated  from  the  maximum  distortion  energy  failure  cri¬ 
terion  using  the  following  expression: 


i 


FIGURE  A. I  Two  Member  FrssJi  and  Typical  Member  Cross-Section 


FIGURE  A. 2  Free  Body  Diagram  for  Two  Member  Frame 
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“ sja 2  *  3-t2,  (A-l) 

where  a  is  the  beading  stress  and  t  is  th«  shear  stress  which  results 
from  twisting  the  members.  Twisting  moments  are  applied  only  at  the 
ends j  therefore,  the  shear  stress  due  to  torsion  is  constant  along 
the  members.  The  bending  moment  varies  along  a  member  but  will  be  a 
tnAxImim  either  at  the  ends  or  under  the  applied  load,  since  the  cross- 
section  is  constant. 

The  failure  stress  will  be  a  maximum  when  the  bending  stress  is 
a  maximum  since  the  torsional  shear  stress  is  constant  along  a  given 
member.  Consequently,  the  maximum  value  of  of  must  occur  at  one  of 
five  points  on  the  two  member  frame.  The  five  possible  points  are  the 
ends  of  each  member  and  wader  the  applied  load.  A  free  body  diagram 
of  member  (l,2)  is  shown  on  Figure  A. 3.  The  twisting  moment  at  any 
section  is 

T  -  Tx  (A-2) 

and  the  bending  moments  at  the  three  critical  points  are: 
at  point  1  (x  *  0) 

M  -  -Mlt  (A-3) 

at  point  4  (x  *  C^) 

M  »  +  R-,0^, 

and  at  point  2  (x  ■  C^g) 

M  »  -Mj.  +  Rici2  “  P^C12  “  ClU^* 

Similarly,  for  member  (2,3)  (see  Figure  A.5),  the  twisting  moment  at 


f 


any  cross-section  i3 


T  »  T, 


(A-U) 


and  the  bending  moments  at  the  two  critical  points  are: 
at  point  3  (x  *  0) 

M  * 

and  at  point  2  (x  *  C^) 

M  .  -M3  +  r3c23. 

The  shear  stress  due  o  transverse  loads  has  been  neglected. 
The  bending  stresses  are  computed  from  the  equation 


(A-5) 


(A-6) 


I 


(A-7) 


where  M  is  the  bending  moment,  h  is  the  height  of  the  rectangular 
section,  and  I  is  the  area  moment  of  inertia  of  the  section.  Using 
the  theory  of  torsion  of  thin  tubes,  the  torsional,  shear  stress  is 
obtained  as 


T 

2At  * 


(A-8) 


where  T  is  the  twisting  moment,  A  is  the  cross-sectional  area,  and  t 
is  the  thickness  of  the  tube.  The  stresses  can  be  calculated  using 
equations  (A-7)  and  (A-6)  if  the  end  reactions  of  the  members  are 


known. 


The  frame  member  shown  in  Figure  A.>  is  subject  to  bending  and 
twisting.  The  strain  energy  due  to  bending  may  be  written  in  terms 
of  the  displacements  of  the  ends  of  the  member. 
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IT 


vhere 


4  * 


V-ortT 

9Z  "1 


U  a  K[62  +  8^2  +  e|  -  3^(8^  +  Sg)  «*  34>23  +  cosset  (A«9) 
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and  X  » 


The  strain  energy  due  to  torsion  is 


,,  GJ  o 
u  a  a  , 


(A*10) 


where  a  is  the  relative  angle  of  tvist  between  the  ends  of  the  Kestber* 
The  potential  energy  of  the  external  load  is 


£1  « 


CL 

-  \  pydx, 
JO 


(A-ll) 


where  both  p  and  y  axe  functions  of  x,  For  the  problems  discussed 
herein  p  is  e  concentrated  load*  The  deflection,  y,  any  be  written 
as  &  function  of  the  end  displactsaento  as  follows: 

x2  -3 

y  *  yl  +  81X  ~  ^281  *  a2  "  TT  *  *6i  ♦  e2  *  2**  ~  *  con8t« 

,  (A-12) 

/ 

Equations  (A-9),  (A~1Q) ,  and  (A-ll)  sre  used  to  c-.btain  the  total  poten¬ 
tial  energy  of  the  fraae  and  the  applied  loads, 

V  °  UT  +  ftj,  (A-13) 

where  U^,  is  the  total  strain  energy,  and  0^  is  the  total  potential 
energy  of  the  external  loads.  By  the  principle  of  stationary  potential 
energy,  the  required  conditions  for  equilibria  are  that 


3V 

3<L 


1  *1,  2,  •••  n. 


(A-l<4) 


where  the  are  generalized  coordinates  (the  three  displaeesent 


i 


components  at  point  2).  Equation  (A-lb)  represents  a  set  of  n  simul¬ 
taneous  lineal'  equations  in  the  n  ucknovn  generalised  coordinates. 

Once  the  deflections  hare  been  determined*  the  end  reactions  may 
be  computed  using  the  following : 


^  *  -K(2ei  ♦  e2  -  30)  + 


(A-15) 


pjj2v 

M0  *  -K(*.  +  20-  -  34)  -  Zz-Z 
212  L2 

SX  "  -  f  (91  +  62  *  S*>  + 


(A-16) 


(A-17) 


and  S2  «  -  (8i  +  82  -  20}  -  — *- 


(A-a8) 


It  should  be  noted  that  it  is  not  necessary  to  write  similar  expres¬ 
sions  for  xhe  twisting  moments.  Due  to  the  relationship  of  action 
and  reaction,  the  twisting  moment  in  one  member  is  just  a  bending 
moment  in  an  adjacent  member. 


Two Member  Frames 


The  generalised  coordinates  are  assigned  as  follows.  The  verti¬ 
cal  displacement  of  point  2  is  q^,  the  joint  rotation  about  the  axis 
of  member  (2,3)  is  q^»  a nd  the  Joint  rotation  about  the  axis  of  member 

C 

(1,2)  is  q^* 

Using  equation  (A-ll)  the  3train  energy  due  to  bending  is 


» ■  14  -  3  ^  H 

C12  *  C12  2  C12 


(A-19) 
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Using  equation  (A-2)  the  strain  energy  due  to  twisting  may  be  written 
as 

®1  0^1  o 

(A-20) 


G12J12  2 
U  *  20^  V 


For  member  (2,3)  the  strain  energy  due  to  bending  is 
U 


•»  4-1  4-i  5 

— g^-~-  Cq2  *  3(~)q,  +  3(~)  3, 

C23  3  °23  3  C23 


(A-21) 


and  the  strain  energy  due  to  twisting  is 

u 

“23  V 


(A-22) 


Equations  (A-ll)  and  (A-12)  are  used  to  obtain  the  potential  energy 
of  the  applied  loads. 


3Cf.  2C3,  C3.  C\ 

a  «  -p(_J -1  „  _»)-  .  P(-i± .  Jdw 
C2  f>3  1  r2  n  ^ 

^12  U12  C-,rt 


(A-23) 

'12  "12 

Adding  the  results  of  equations  (A- 19)  through  (A-23),  the  total 
potential  energy  of  the  frame  is 
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(A-21*) 
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When  the  principle  of  stationary  potential  energy  is  applied,  three 
equations  in  q^  q^t  and  q3  are  obtained  from  equation  (A-24 ) 


|i-  -  [12E12I12  ,  6£12X12  ^  _  ffg£ 
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(A-25) 
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Equations  (A-7)  through  (A-10)  are  nov  used  to  obtain  the  end 
loads  on  the  ambers.  Only  the  loads  at  points  1  and  3  (see  Figure 
A-2)  are  required  to  determine  the  stresses. 


%  -  -  ■  3  gSa  *  c^)2 

°12  *  <%> 


(A-29) 
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•Tiie  stresses  are  eoaputed  using  equations  (A— 7 )  and  (A-8) .  As 
can  be  seen  from  these  equations,  and  equations  (A-2)  through  (A-6), 
the  stresses  will  he  linear  functions  of  the  end  reactions.  Expressed 

in  matrix  fora  the  stresses  say  he  eoaputed  from  the  following  equa¬ 
tion: 

S  »  CM  +  P,  (A— >36 ) 

vhere  M  is  the  vector  of  end  reactions  defined  previously.  The  matrix 
C  and  the  vectors  P  and  S  are  defined  below. 
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stresses  on  member  (2*3)  at  points  2  said  3.  and  are  the  shear¬ 
ing  stresses  due  to  torsion  in  member  (1*2)  and  (2*3). 

A. 3  Application  to  Three  Member  Frames 

Equations  (A-l)  through  (A-18)  may  be  applied  to  the  three  member 
frame  by  following  the  procedures  outlined  in  Section  A. 2.  The 
resulting  matrix  equations*  (A-28),  (A-35),  and  (A-36),  are  repeated 
here. 


A<i  »  P 

(A-28) 

M  ■  BQ  +  F  • 

(A-35) 

S  «  CM  +  P 

(A-36) 

These  equations  may  be  used  directly  by  defining  the  necessary  vectors 

and  matrices  l’or  the  three  member  frame  shown  in  Figure  A. 6, 

The  deflections,  Q*  are  defined  as  follows: 
qj ,  is  the  vertical  displacement  of  Joint  3; 
qg,  the  rotation  of  Joint  3  about  axis  (3,4); 

q^»  the  rotation  of  Joint  3  about  axis  (1,3); 

q^,  the  vertical  displacement  of  Joint  5; 
q^*  the  rotation  of  Joint  5  about  axis  (3*5); 

and  <j£»  the  rotation  of  Joint  5  about  axis  (5*7). 

The  vectors  M  and  S  are  defined  using  the  free  body  diagram  in 
Figure  A, 7. 
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FIGURE  A. 6  Three  Member  Frame 
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FIGURE  A.  7  Free  Body  Diagram  of  the 
Three  Member  Frame 
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1*7 


and 


where 


is  the  heeding  stress  at  point  1  on  member  (l,3)t 
Og,  the  bending  stress  at  point  2  in  member  (1,3), 

a^,  the  bending  stress  point  3  in  member  (1,3); 

c^,  the  torsional  stress  in  member  (1,3); 
act  the  bending  stress  at  point  7  in  membej  (5,7); 

o^,  the  bending  stress  at  point  6  in  member  (5,'T); 

o^.,  the  bending  stress  at  point  5  in  member  (5,7); 

Oq ,  the  torsional  stress  in  member  (5,7); 


DOCUMENT  CONTROL  DATA  -  R  &  D  • 

_ c*  irjtxlnf  annotation  must  bo  mnterid  when  tho  owrmii  rcpott  lx  ci*t»l(lmd) 

1  OR'CiNATlNO  ACTIVITY  (Corpoff  0Uthor)  j2*.  REPORT  SECURITY  C  LASS1  FIG  ATION 

The  University  of  Iowa  [ujggssified - . 


|3.  REPOR1  TITLE 


OPTIMUM  DESIGN  OF  SPATIAL  STRUCTURES 


VS  *•  DESCRIPTIVE  notes  (Typs  ol  report  and  Inclusive  detee) 

N  Technical  Report _ 

s-  AOTHOPW)  (First  netne,  middls  Inltlel,  lest  neme) 

n5"  Donald  L.  Bartel 

A _ 

6-  REPORT  DATE 

August  1969 _ 

ML  CONTRACT  OR  GRANT  NO. 

DAAF 03-69-C-0014 

be  PROJECT  NO. 


17A.  TOTAL  NO.  OP  PAGES 


|7 b.  NO.  OF  REFS 


122 _ I  27 

9a.  ORIGINATOR* S  REPORT  NUkPSER(S) 

Project  Themis  Report  No.  12 


|  10.  DISTRIBUTION  STATEMENT 


Oh.  OTHER  REPORT  NO(f>  (Any  oihmr  number*  thmt  mmy  6«  ««»/0n«cf 
thin  report) 


Distribution  of  this  report  is  unlimited 


It.  SUi  PLEMENTARY  NOTES 


t2.  SPONSORING  MILITARY  ACTIVITY 

U.S.  Army  Weapons  Command 
R*E  Directorate 


M3*  ABSTRACT 


This  research  presents  a  systematic  approach  to  the  optimal  design  of 
spatial  structures  for  minimum  weight  subject  to  constraints  on  stress  and  geometry. 
The  optimization  procedures  discussed  are  general,  and  may  be  applied  to  structures 
which  can  be  analyzed  by  matrix  displacement  or  finite  element  methods. 

Two  methods  of  mathematical  programming  are  applied  to  obtain  a  minimum 
weight  design.  The  first  is  the  sequential  unconstrained  minimization  technique 
(SUMT)t  and  the  second  is  the  method  of  constrained  steepest  descent  with  state 
equations  (CSDS).  Both  of  these  techniques  require  derivatives  of  the  objective 
and  constraint  f unct ions  to  improve  estimates  of  the  optimum  design.  In  many 
structural  problems,  it  is  very  difficult  or  impossible  to  compute  these  deriv¬ 
atives  exactly;  existing  structural  analysis  algorithms  are  generally  not 
equipped  to  compute  these  derivatives.  In  order  co  take  full  advantage  of 
existing  analysis  capability,  the  programming  techniques  in  this  research  h3ve 
been  developed  assuming  that  sudi  derivatives  are  not  available. 

Optimal  structural  design  problems  arc  characterized  by  an  objective  function 
(the  weight),  state  variables  (the  stresses  and  deflections),  design  variables, 
state  equations  (the  structural  analysis),  and  constraints  which  may  be  functions 
of  the  design  and  state  variables.  When  the  state  eq,T^tions  are  used  to  write 
all  of  the  constraints  as  functions  of  the  design  variables,  a  nonlinear  pro¬ 
gramming  problem  results.  The  sequential  unconstrained  minimizagion  technique 

reduces  the  constrained  nonlinear  programming  problem  to  a  semiPnn. 
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problems  which  can  be  solved  using  existing  unconstrained  minimization 
techniques,  A  SUMT  program  was  written  for  this  research  using  Powell's 
method  of  unconstrained  minimization  without  derivatives.  The  required 
minimization  of  a  function  along  a  line  uses  a  combination  of  a  Fibonacci 
search  (to  bracket  the  minimum)  and  a  quadratic  approximation  of  the 
minimum. 

The  method  of  constrained  steepest  descent  diffsrr  from  the  usual 

nonlinear  programming  problem  in  that  the  state  equations  and  the  state 

variable  constraints  appear  explicitly  in  the  formulation.  This  provides 

a  natural  matching  of  the  essential  features  of  the  design  problem  and  the 

method  used  to  obtain  its  solution.  The  design  problem  is  linearised 

about  a  candidate  design  and  the  desired  improvement  in  the  design  variables, 

T  -1  2 

5x,  is  required  tc  be  small  by  demanding  that  6x  w  6x  =  £  r,  where  £  is  a 
small  number  and  w  is  a  positive  definite  weighting  matrix.  The  Kuhn-Tucker 
necessary  conditions  are  then  applied  to  the  resulting  nonlinear  problem. 

As  a  direct  consequence,  ox  is  specified  in  terms  of  two  components;  6x^ 
which  reduces  the  objective  function  consistent  with  the  constraints,  and 
6x2  which  directs  the  search  for  a  minimum  back  to  the  feasible  region  if 
constraints  have  been  violated.  The  method  was  applied  using  both  exact 
and  approximate  derivatives,  so  that  its  effectiveness  when  derivatives 
are  not  available  could  be  assessed, 

A  spatial  structures  which  occurs  frequently  in  practice  is  the  plane  frame 
with  out-of-plane  leads.  Although  such  structures  are  generally  made  up  of 
relatively  few  members,  they  may  have  many  design  variable.;  since  several  design 
parameters  must  be  specified  for  each  member.  The  programming  methods  were  applied 
to  a  number  of  two  and  three  member  frames  of  this  type.  From  the  results,  it 
appears  that  CSDS  has  significant  advantages  over  SUMT  both  in  terms  of  compu¬ 
tational  time  and  the  number  of  times  that  candidate  designs  must  be  analyzed. 

The  results  also  show  that  CSDS  performs  as  well  when  derivatives  are  approx¬ 
imated  as  it  does  when  they  can  be  computed  exactly.  The  effectiveness  of  SUMT 
is  reduced  significantly  if  the  derivatives 
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